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Abstract 


The  main  goal  of  this  project  is  to  estimate  theoretically  the  r  "trail  or  effective  constitutive 
properties  of  nonlinear  composite  materials  undergoing  large  deformations.  Two  types  of 
large  deformations  are  of  interest:  large  elastic  deformations,  and  large  viscous 
deformations.  The  proposed  method  is  to  apply  variational  principles  that  are  under 
development  to  characterize  the  range  of  the  effective  properties  given  partial  statistical 
information  about  the  microstructure  (such  as  the  volume  fractions  of  the  phases).  For 
some  particular  microstructures  of  interest  exact  estimates  may  be  given.  Significant 
progress  was  made  over  the  first  year  with  the  development  of  a  new  variational  principle 
allowing  the  estimation  of  the  effective  properties  of  a  given  nonlinear  composite  in  terms 
of  the  etfective  properties  of  linear  composites  (which  are  assumed  to  be  known).  The 
potential  significance  of  this  work  derives  from  its  simplicity  allowing  the  application  of  a 
large  body  of  prior  research  on  linear  composites  to  nonlinear  composites.  This  method  has 
been  applied  to  the  case  of  large  viscous  deformations,  and  some  results  for  particular 
materials  systems  have  already  been  reported  in  the  pertinent  literature. 


Research  goals 

The  main  goal  of  this  project  is  to  estimate  the  overall  or  effective  constitutive  properties 
of  nonlinear  composite  materials  undergoing  large  deformations.  Two  types  of  large 
deformations  are  of  particular  interest:  large  elastic  deformations,  corresponding  to 
materials  such  as  polymeric  composites,  rubber  foams  and  solid  rocket  fuel  composites; 
and  large  viscous  deformations,  corresponding  to  the  high-temperature  creeping,  or  to  the 
dynamic  plastic  deformation  of  metals. 

Background 

The  first  and  only  available  rigorous  procedure  thus  far  for  estimating  the  overall 
constitutive  properties  of  nonlinear  random  composites  was  given  by  Talbot  &  Willis 
(1985).  This  procedure  is  based  on  an  extension  of  the  well-known  Hashin-Shtrikman 
(HS)  variational  principles  to  a  class  of  nonlinear  materials.  Ponte  Castaneda  &  Willis 
(1988)  applied  this  procedure  to  nonlinearly  viscous  materials,  and  gave  the  first  rigorous 


bounds  and  self-consistent  (SC)  estimates  for  the  effective  properties  of  composites  in  a - 

power-law  creep.  Ponte  Castaneda  (1989)  extended  the  work  of  Talbot  &  Willis  to  finite 
elasticity,  and  used  the  resulting  theory  to  provide  the  first  bounds  and  estimates  for  ;,/ 

nonlinearly  elastic  composites.  ty  Codes 

e^:d/or 

Special 


□  □ 


Accomplishments 

Introduction 

The  research  carried  out  under  the  sponsorship  of  AFOSR  in  the  seventeen  months  has 
brought  about  an  important  development,  reported  in  Ponte  Castaneda  (1990a  &  b).  In 
these  works,  a  new  variational  structure  is  introduced  that  yields  the  following  important 
result  for  the  effective  properties  of  isotropic  nonlinear  composites:  given  the  overall 
properties  of  a  "comparison"  linear  composite  with  the  same  microstructural 
distribution  of  phases  as  the  nonlinear  composite,  it  is  possible  to  give  an  estimate  (in 
the  form  of  a  bound)  for  the  overall  properties  of  the  nonlinear  isotropic  composite. 
Thus,  available  bounds  and  estimates  for  linear  composites  can  be  "translated"  into 
corresponding  bounds  and  estimates  for  nonlinear  composites.  The  new  structure  has  the 
following  advantages:  it  is  more  general  than  the  corresponding  structure  of  Talbot  & 
Willis,  because  it  not  only  yields  HS  bounds  and  SC  estimates,  but  it  generalizes  to  the 
nonlinear  case  any  other  type  of  bound  or  estimate  that  may  be  available  for  the  linear 
composite;  it  is  easier  to  implement,  because  it  only  requires  the  calculation  of  the  extreme 
values  of  two  appropriately  defined  functions;  and,  finally,  it  gives  in  some  cases  stronger 
(in  other,  identical)  results  than  the  method  of  Talbot  &  Willis. 

The  method  has  been  applied  to  three  different  material  systems.  Porous  and  rigidly 
reinforced  materials  with  pure-power  law  and  Ramberg-Osgood  matrix  behaviors, 
respectively,  are  considered  in  Ponte  Castaneda  (1990a)  and  Ponte  Castaneda  and  de 
Botton  (1990).  Alternatively,  Ponte  Castaneda  (1990b)  considered  a  family  of  brittle/ductile 
composites  containing  a  linear  and  a  nonlinear  phase.  In  addition  to  low-temperature 
plasticity  (in  its  deformation  theory  form),  these  results  can  also  be  re-interpreted  in  the 
context  of  high-temperature  creep  and  high  strain-rate  viscoplasticity. 

The  method 

Consider  an  n-phase  composite  occupying  a  domain  Q  (normalized  to  have  unit  volume), 
with  each  phase  occupying  a  sub-domain  (r  =  1,2,...,  n),  and  let  the  stress  potential, 
f/(o,x),  be  expressed  in  terms  of  the  n  homogeneous  phase  potentials,  Uir\a),  via 

U(a,x)=XzM(x)U(r)(a), 

r=l 


(1) 


where 


Jl  ifxea(r) 

[O  otherwise 


(2) 


is  the  indicator  function  of  phase  r.  The  phases  are  assumed  to  be  isotropic,  so  that  the 
potentials  Uir)(a)  depend  on  the  stress  a  only  through  its  three  principal  invariants.  Here, 
we  will  further  assume  that  the  dependence  is  only  through  two  of  the  invariants,  namely. 


the  mean  stress,  a.n  =  -trc,  and  the  effective  stress,  <r,  =  J— S-S,  where  S  is  the 


deviator  of  O. 


The  effective  behavior  of  the  composite  material  is  defined  in  terms  of  the  effective  energy, 
0(a),  that  arises  due  to  the  uniform  traction  boundary  condition 

aijnj  =  aiJn),  xedQ,  (3) 

where  dC2  denotes  the  boundary  of  the  composite,  n  is  its  unit  outward  normal,  and  a  is  a 
given  constant  symmetric  tensor.  The  average  stress  in  the  composite  is  then  precisely  a 
and  it  follows  from  the  principle  of  minimum  complementary  energy  that 


where 


0(a)  =  min  U(a), 

aeS(a\ 


(4) 


t/(o)=[  U(a,x)dV 
Jn 

is  the  complementary  energy  functional  of  the  problem,  and 

S(a)  =  {al Gjjj  =  0  in  Q,  and  =  oijnl  on  dQ) 
is  the  set  of  statically  admissible  stresses.  Thus,  if  E  denotes  the  average  strain  over  the 
composite,  it  can  be  readily  shown  that 


e  -—-(a),  (5) 

y  dau 

which  yields  an  effective  constitutive  relation  for  the  composite  in  terms  of  the  average 
variables  a  and  E.  Given  this  connection  between  the  effective  potential  for  the  composite 
U(a)  and  the  effective  stress/strain  relation,  it  makes  sense  to  seek  information  on  0(a). 


Next  we  make  use  of  a  linear  heterogeneous  "comparison"  material,  with  effective 
properties  that  can  be  characterized  in  terms  of  bounds  and  estimates,  to  obtain 
corresponding  bounds  and  estimates  for  the  effective  properties  of  a  nonlinear  composite. 


To  this  end,  we  introduce  the  quadratic  potential 


U(a,x)  =  f  z<r)(x)U<r,(o)  =  ~J—  a2  + 
%  6/i(x) 


_ 1 

2k{x) 


oi. 


(6) 


such  that /i(x)  =  ^^(,')(x)/i<'',>0,  and  k(x)  =  ^%lr)(x)Kir)  >  0,  with  the  /i(r)  and  k(r) 

r= 1  r=l 

constant,  corresponding  to  a  linear  isotropic  composite  with  the  same  phase  distribution 
(the  same  indicator  functions)  as  the  nonlinear  composite. 

Then,  if  the  nonlinearity  in  the  potential  of  the  original  composite  is  stronger  than  quadratic 
as  the  norm  of  the  stress  becomes  large,  it  makes  sense  to  define  the  set  of  functions 

V(r\ju(r),  k(r))  =  mzx[0(r\o)  -  t/(r)(CT)}.  (7) 

O 

It  follows  that, 


where 


0(a)  >  max  {< U(o)-V(fi,k )). 


U(a)=  min  U(a) 

oeS(3) 

is  the  effective  potential  of  the  linear  composite,  and 


(8) 

(9) 


V(fi  ,ic)  =  ^cir)Vlr)(fi(r\  klr)),  (10) 

r=l 

is  expressed  in  terms  of  the  volume  fractions  of  each  phase, 

c(r)  =  J 'nxir\x)dV. 

The  details  of  the  derivation  of  this  result  are  given  in  Appendix  A  (Ponte  Castaneda  1990a) 
in  pages  6  through  9.  We  note  however  that  expression  (8)  allows  the  estimation  of  the 
effective  properties  of  the  nonlinear  composite  in  terms  of  the  effective  properties  of  a 
family  of  linear  composites  with  elastic  moduli  fi{r)  and  k{r). 

Usually,  however,  it  is  not  possible  to  find  0(o)  explicitly,  but  instead  bounds  and 
estimates  may  be  available  for  0(o).  If  we  have  a  lower  bound  (such  as  a  Hashin- 
Shtrikman  lower  bound)  U_(o)  for  U(o),  such  that 

U(o)>U_(o),  (11) 

then,  replacing  U(o)  by  0_(o)  in  equation  (8)  for  0_(a)  yields  a  lower  bound  for  0(a). 


On  the  other  hand,  if  we  only  have  an  estimate  (such  as  a  self-consistent  estimate)  Ut(o) 
for  £/(o),  then 

0  (a)  =  max  [Ut(a)~  V(fi,k)}  (12) 

/iw.rw>  0 

would  provide  only  an  estimate  for  0(a). 

We  note  that  the  prescriptions  (8)  and  (12)  lead  to  convex  expressions  for  the  bounds  and 
estimates  of  the  effective  potential  of  the  nonlinear  composite,  provided  that  the 
corresponding  bounds  and  estimates  for  the  linear  composite  are  convex,  which  is  in  turn 
guaranteed  assuming  that  (i  and  k  >  0.  This  is  a  desirable  feature,  because  the  effective 
potential  of  the  composite  is  known  to  be  convex. 

An  alternative  derivation,  and  a  stronger  version,  of  this  variational  procedure  is  given  in 
Appendix  C  (Ponte  Castaneda  1990b)  in  pages  4  to  6.  The  advantage  of  the  alternative 
derivation  is  that  it  allows  the  characterization  of  the  circumstances  under  which  the 
inequality  in  equation  (8)  turns  into  equality.  This  idea  is  important  in  the  context  of 
assessing  the  strength  of  the  bounds. 

Results 

The  general  procedure  was  applied  in  Ponte  Castaneda  (1990a)  to  an  isotropic  porous 
material  with  incompressible  behavior  for  the  matrix  with  potential  described  by 

U('\a)  =  f(o,).  (13) 

The  result  of  this  calculation  are  a  Hashin-Shtrikman  (H-S)  lower  bound  for  U(a),  given 
by 

0.(0)  =  c{1)/(s),  (14a) 

with 

5  =  -jjy  ^(1  +  % c(2)) a]  +  %  c(2)  a2n ,  (14b) 

and  a  self-consistent  (S-C)  estimate  given  by 

0,(a)  =  c(1)/(i), 


(15a) 


with 


(15b) 


The  details  of  the  derivation  of  these  results  is  given  in  pages  9  to  12  of  Appendix  A. 
Results  are  also  given  in  this  reference  for  a  rigidly  reinforced  material  with  matrix  behavior 
characterized  by  relation  (13),  and  for  a  general  two-phase  incompressible  composite.  The 
above  results  are  also  specialized  in  Ponte  Castaneda  (1990a)  to  pure  power-law  behavior 
for /.  Results  for  appropriately  defined  low-triaxiality  "shear”  and  high-triaxiality  "bulk" 
moduli  for  the  porous  material  are  given  in  Figures  2  and  3  of  Appendix  A  for 
representative  values  of  the  hardening  parameter  n  as  functions  of  the  porosity.  Analogous 
results  are  given  in  Figures  3  and  4  for  the  effective  "shear"  modulus  of  the  two-phase 
incompressible  composite  and  the  rigidly  reinforced  material,  respectively.  The  new  results 
are  compared  with  previously  available  results  such  as  dilute  and  self-consistent  estimates. 
For  the  porous  material,  the  bounds  are  found  to  be  the  best  bounds  available  (better  than 
the  corresponding  bounds  of  Ponte  Castaneda  and  Willis,  1988  obtained  using  the 
nonlinear  extension  of  the  Hashin-Shtrikman  variational  principle  of  Talbot  and  Willis, 
1984).  In  addition  to  their  intrinsic  value,  these  bounds  are  also  essential  in  characterizing 
the  range  of  validity  of  numerical  calculations  based  on  the  dilute  approximation.  Thus,  it  is 
found  (see  Figures  2)  that  the  dilute  calculations  of  Duva  and  Hutchinson  (1984)  for  the 
low-triaxiality  modulus  are  reasonably  good  for  porosities  up  to  40%.  On  the  other  hand, 
their  corresponding  calculations  for  the  high-triaxiality  bulk  modulus  (Figures  3)  have  a 
very  small  range  of  validity  (less  than  1%  for  practically  important  values  of  the  hardening 
parameter).  In  Figure  4,  the  new  self-consistent  estimates  for  the  rigidly  reinforced  material 
compare  favorably  with  the  corresponding  differential  self-consistent  estimates  of  Duva 
(1984). 

In  Ponte  Castaneda  and  de  Botton  (1990),  the  above  procedure  is  specialized  to  a  porous 
material  with  "linear  plus  power  hardening"  (Ramberg-Osgood)  behavior  for  the  nonlinear 
phase.  Results  are  given  for  the  effective  potential  of  the  composite  as  a  function  of  the 
average  stress  and  for  representative  stress/strain  relations  for  the  porous  material  in 
Figures  1  and  2  of  Appendix  B  for  different  values  of  the  hardening  parameter  and  the 
porosity.  The  results  are  compared  with  results  by  Willis  (1990),  and  the  new  results  are 
found  to  be  superior.  In  particular,  the  physically  unrealistic  discontinuous  behavior  of 
Willis'  model  between  the  linear  and  nonlinear  domains  at  the  larger  values  of  the  triaxiality 


(see  the  continuous  lines  in  Figure  Id)  is  not  observed  in  the  new  results  (there  is  a  smooth 
transition  between  the  linear  and  plastic  domains  for  the  dashed  lines). 

Finally,  in  Ponte  Castaneda  (1990b)  results  in  the  foim  of  Hashin-Shtrikman  bounds  are 
given  for  composites  containing  linear  and  nonlinear  (c.f.  equation  (13))  phases.  As  it  turns 
out,  some  of  these  bounds  can  be  shown  to  be  optimal  (best  possible  bounds  given  the 
volume  fractions  of  each  phase)  in  some  cases.  This  is  accomplished  by  identifying  special 
microstructures  (called  sequentially  laminated  composites)  attaining  the  bounds.  Results 
specialized  to  pure  power-law  behavior  for  the  nonlinear  material  are  given  in  Figures  2  and 
3  of  Appendix  C.  Figures  2  give  bounds  for  an  appropriately  normalized  form  of  the 
effective  energy  of  the  composite  versus  an  appropriately  normalized  measure  of  the 
average  stress  for  representative  values  of  the  hardening  parameter  n  at  small  enough  stress 
levels.  Figures  3  give  the  corresponding  results  for  larger  values  of  the  average  stress.  In 
Figures  4,  some  results  are  given  for  a  linear  plus  power  (Ramberg-Osgood)  behavior  for 
the  nonlinear  phase.  These  results  for  the  bounds  are  the  first  of  its  kind,  and  thus  we  have 
not  been  able  to  compare  them  with  other  results. 

Research  plans  for  the  future 

The  work  proposed  initially  in  this  project  has  been  briefly  interrupted  by  the  move  of  the 
PI  and  the  graduate  students  under  his  supervision  from  the  Johns  Hopkins  University  to 
the  University  of  Pennsylvania.  This  move  is  expected  to  enhance  the  research  activities  of 
the  PI  through  his  interactions  with  the  NSF-sponsored  Materials  Research  Laboratory  at 
Penn.  The  PI  has  requested  continuation  of  the  project  at  Penn,  and  this  has  been 
recommended  by  AFOSR  for  sponsorship  starting  January  1,  1991.  Encouraged  by  the 
results  obtained  at  Hopkins,  we  plan  to  continue  the  work  proposed  in  the  original 
proposal.  First,  we  would  like  to  address  anisotropic  composites  (such  as  fiber-reinforced 
materials),  and  then  later  nonlinearly  elastic  (with  finite  deformations)  composites.  The 
extension  to  anisotropic  composites  is  already  in  progress,  with  the  assistance  of  one  of  my 
graduate  students,  Mr.  Gal  de  Botton,  who  is  being  supported  by  the  Mechanical 
Engineering  Department  in  the  interim  period  until  the  new  project  becomes  active.  This 
work  is  expected  to  yield  important  new  results,  and  it  will  serve  to  introduce  Mr.  de 
Botton  to  this  area  of  research.  Additionally,  another  graduate  student,  on  a  parallel  project 
sponsored  by  NSF,  will  attempt  to  concentrate  on  the  special  case  of  porous  materials 
(composites  with  a  vacuous  phase).  Finally,  the  non-trivial  extension  to  nonlinearly  elastic 
composites  is  expected  to  occupy  our  efforts  for  the  most  part  of  the  next  nineteen  months. 
The  complications  in  this  last  problem  are  related  to  the  strong  nonlinearities  that  arise  due 
to  the  large  deformation  kinematics,  but  it  is  expected  that  the  ideas  developed  in  Ponte 


Castaneda  (1989)  will  also  be  useful  in  this  work,  leading  to  improved  results  for  the 
effective  properties  of  such  composites. 
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Abstract 

A  new  variational  structure  is  proposed  that  yields  a  prescription  for  the  effective  energy  potentials  of 
nonlinear  composites  in  terms  of  the  corresponding  energy  potentials  for  linear  composites  with  the  same 
microstructural  distributions.  The  prescription  can  be  used  to  obtain  bounds  and  estimates  for  the  effective 
mechanical  properties  of  nonlinear  composites  from  any  bounds  and  estimates  that  may  be  available  for 
the  effective  properties  of  linear  composites.  The  main  advantages  of  the  procedure  are  the  simplicity  of 
its  implementation,  the  generality  of  its  applications  and  the  strength  of  its  results.  The  general  prescription 
is  applied  to  three  special  nonlinear  composites :  a  porous  material,  a  two-phase  incompressible  composite 
and  a  rigidly  reinforced  material.  The  results  are  compared  with  previously  available  results  for  the  special 
case  of  power-law  constitutive  behavior. 


1 .  Introduction 

The  prediction  of  the  effective ,  or  overall,  constitutive  behavior  of  composite  solid 
materials  is  both  a  practically  and  theoretically  important  problem,  which  draws  input 
from  many  different  disciplines,  including  material^  science,  mechanics  and 
mathematics.  This  paper  deals  with  the  theoretical  prediction  of  the  effective  mechan¬ 
ical  properties  of  heterogeneous  materials  with  nonlinear  phases  that  are  perfectly 
bonded  to  each  other,  and  isotropic.  To  make  sense  of  the  notion  of  effective  properties 
for  the  composite,  the  size  of  the  typical  heterogeneity  is  generally  assumed  to  be 
small  compared  to  the  size  of  the  specimen  and  the  scale  of  variation  of  the  applied 
loads.  It  is  further  assumed  that  the  effect  of  the  interfaces  is  negligible,  so  that  the 
effective  properties  of  the  composite  are  essentially  derived  from  the  bulk  behavior  of 
the  constituent  phases. 

The  corresponding  theory  for  linear  composites  is  fairly  well  developed,  including 
different  approaches  to  the  problem  with  varying  degrees  of  mathematical  sophis¬ 
tication  and  physical  relevance.  Thus,  exact  estimates  have  been  determined  for  the 
effective  properties  of  ad  hoc  models  of  composites ;  rigorous  variational  bounds  have 
been  given  for  the  properties  of  random  composites;  and  precise  definitions  and 
explicit  “homogenization”  formulae  have  been  proposed  for  the  properties  of  periodic 
composites.  Appropriate,  but  by  no  means  exhaustive,  references  dealing  with  the 
linear  theories  are  provided  by  the  review  articles  of  Willis  (1982.  1983)  and  Kohn 
(1989),  and  by  the  monographs  of  Christensen  (1979)  and  Sanchez-Palencia 
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(1980).  By  contrast,  in  the  context  of  nonlinear  composites,  most  of  the  results  are 
based  on  ad  hoc  models,  such  as  dilute  and  self-consistent  models.  For  instance,  Duva 
and  Hutchinson  (1984)  (referred  to  as  DH  fo-thfr-seewSh  gave  dilute  estimates  for 
the  effective  properties  of  a  nonlinear  porous  material,  and  Duva  (1984)  proposed 
self-consistent  estimates  of  the  differential  type  for  a  rigidly  reinforced,  nonlinear 
material.  To  the  knowledge  of  the  author,  the  first  and  only  contribution  so  far  dealing 
with  the  calculation  of  rigorous  bounds  for  the  effective  properties  of  nonlinear 
composites  is  provided  by  the  work  of  Talbot  and  Willis  (1985)  (referred  to  as  TW), 
which  was  introduced  in  Willis  (1983).  These  authors  extended  the  well-known 
Hashin-Shtrikman  variational  principles  to  include  nonlinear  constitutive  behavior, 
and  their  methods  have  been  applied  to  a  number  of  examples  in  different  physical 
contexts.  For  example,  Ponte  Castaneda  and  Willis  (1988)  (referred  to  as  PCW), 
and  more  recently  Willis  (1989),  have  determined  bounds  and  estimates  for  the 
effective  properties  of  nonlinearly  viscous  (or  infinitesimally  elastic)  materials.  Also, 
Ponte  Castaneda  (1989)  jgam  provided  extensions  of  the  minimum  complementary 
energy  and  the  Talbot-Willis  variational  principles  to  finite  elasticity,  that  allowed 
the  calculation  of  bounds  and  estimates  for  the  effective  properties  of  a  broad  class  of 
nonlinearly  elastic  composites.  Additional  developments  are  provided  by  appropriate 
extensions  of  the  periodic  homogenization  formula  by  Marcellini  (1978)  for  prob¬ 
lems  with  convex  energy  densities,  and  by  Muller  (1987)  for  finite  elasticity  (with  a 
non-convex  energy  density). 

In  this  work,  we  propose  an  alternative  variational  structure  that  allows  the  esti¬ 
mation  of  the  effective  energy  densities  of  nonlinear  composites  in  terms  of  the 
corresponding  information  for  linear  composites  with  the  same  microstructural  dis¬ 
tribution.  Although  the  procedure  has  application  to  problems  in  other  physical 
contexts  that  have  a  variational  characterization  (see  Ponte  Castaneda,  1990,  for  an 
application  in  conductivity),  here  we  will  study  the  specific  application  of  the  theory 
to  composite  materials  with  constitutive  behavior  characterized  by  nonlinear  viscosity, 
or  by  the  mathematically  analogous  theory  of  nonlinear  infinitesimal  elasticity.  Effec¬ 
tive  properties  are  defined  in  Section  2  by  means  of  the  principle  of  minimum  potential 
energy,  and  its  dual  counterpart,  the  principle  of  minimum  complementary  energy. 
The  new  structure  is  developed  in  Section  3,  where  it  is  shown  that  a  general  bound 
(or,  alternatively,  an  estimate)  for  the  linear  composite  can  be  translated  into  a  bound 
(or  estimate)  for  the  nonlinear  composite.  In  Section  4,  the  general  procedure  is 
applied  to  three  particular  cases  of  general  interest:  a  porous  material,  a  two-phase 
incompressible  composite  and  a  material  reinforced  by  rigid  inclusions.  For  each  of 
these  composites,  we  give  bounds  and  estimates  for  their  effective  properties.  In 
Section  5,  the  results  are  specialized  further  to  phases  with  a  power-law  type  of 
constitutive  behavior,  and  the  results  are  discussed  and  compared  with  previously 
available  results.  Finally,  in  Section  6  some  general  conclusions  are  drawn. 


2.  Effective  Properties 

We  are  interested  in  estimating  the  effective,  or  overall,  properties  of  composites 
with  nonlinear  material  behavior.  By  a  "composite"  we  mean  an  idealized  material 
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that  corresponds  to  the  limit  of  a  sequence  of  heterogeneous  materials  with  two 
distinct  length  scales :  one  microscopic  /  corresponding  to  the  size  of  the  heterogeneity, 
and  one  macroscopic  L  corresponding  to  the  size  of  the  specimen  of  interest  and  the 
scale  of  variation  of  the  boundary  conditions.  The  effective  behavior  of  the  composite 
is  then  obtained  by  considering  the  limit  of  the  behavior  of  the  sequence  of  materials 
as  the  ratio  of  scales  s  =  IjL  tends  to  zero.  The  study  of  the  definition  and  existence 
of  such  limit  properties  is  called  homogenization  theory,  and  it  is  an  area  of  current 
interest  in  the  general  mathematics  community  (Kohn,  1989).  However,  for  the 
purposes  of  this  work,  it  will  not  be  necessary  to  introduce  this  formalism ;  we  can 
always  rely  on  our  physical  intuition  to  understand  the  concept  of  effective  properties, 
and  in  the  analysis  that  follows,  it  will  suffice  to  take  our  composite  to  be  a  het¬ 
erogeneous  material  with  very  small,  but  finite  microscale.  The  effective  properties 
are  then  understood  in  an  approximate  sense. 

Consider  an  tf-phase  composite  occupying  a  domain  fi,  with  each  phase  occupying 
a  sub-domain  £2'"  (r  =  1 , 2, ....  n),  and  let  the  stress  potential,  U{<t,  x),  be  expressed 
in  terms  of  the  n  homogeneous  phase  potentials,  U{'\a),  via 

l/(M)=il",(x)f/wW,  (2.1) 

r-  I 


where 


X(,)(x) 


1  ifxefr1 

0  otherwise 


(2.2) 


is  the  characteristic  function  of  phase  r.  The  phases  are  assumed  to  be  isotropic,  so 
that  the  potentials  (7w(<r)  depend  on  the  stress  a  only  through  its  three  principal 
invariants.  Here,  we  will  further  assume  that  the  dependence  is  only  through  two  of 
the  invariants,  namely,  the  mean  stress 

am  -  3  tr  a. 


and  the  effective  stress 


a,  =  yisii, 

where  S  is  the  deviator  of  a. 

Then,  the  strain  tensor  (or  the  strain-rate  tensor,  depending  on  whether  we  are 
dealing  with  nonlinear  infinitesimal  elasticity  or  viscosity)  e,  which  is  required  to 
satisfy  the  compatibility  relations 

=  i(H/., +uu)  (2.3) 

in  terms  of  the  displacement  (or  velocity  field)  u,  is  related  to  the  stress  a ,  satisfying 
the  equilibrium  equations 


via  the  constitutive  relation 


=  0, 


(2.4) 


(2.5) 
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The  commas  in  (2.3)  and  (2.4)  denote  differentiation,  and  the  summation  convention 
has  also  been  used  in  (2.4).  We  assume  that  the  phases  are  perfectly  bonded,  so 
that  the  displacement  (or  velocity)  is  continuous  across  the  interphase  boundaries. 
However,  the  strains  and,  therefore,  the  stresses  may  be  discontinuous  across  such 
boundaries,  and  hence  (2.4)  must  be  interpreted  in  a  weak  sense,  requiring  continuity 
of  the  traction  components  of  the  stress  across  the  interphase  boundaries.  Also,  at 
least  one  of  the  phase  potentials  is  assumed  to  be  non-quadratic  in  the  stress,  so  that 
the  constitutive  response  of  the  material  as  given  by  (2.5)  is  genuinely  nonlinear. 

The  statement  of  the  problem  is  completed  by  the  selection  of  an  appropriate 
boundary  condition : 


Gijnj  ~  Gijnj'  xecQ,  (2.6) 

where  c£2  denotes  the  boundary  of  the  composite,  n  is  its  unit  outward  normal,  and 
6  is  a  given  constant  symmetric  tensor.  This  uniform  constraint  condition  has  some 
useful  properties,  discovered  by  Hill  (1963).  Let 


and 


a  = 


e  = 


1 

1 


<x(x)  dV, 


a(x)dK 


(2.7) 


(2.8) 


denote  the  respective  averages  of  the  actual  stress  and  strain  fields  in  the  composite. 
Here,  the  scale  of  has  been  normalized  to  have  unit  volume.  Then,  we  have  that 

a  =  a,  (2.9) 


and 


(u  g>  n  +  n  ®  u)  dS.  (2.10) 

2  Jen 

This  means  that  the  average  stress  a  in  the  composite  is  precisely  a,  and  that  the 
average  strain  i  can  be  “measured”  in  terms  of  the  boundary  displacements. 

The  third  property  makes  use  of  the  principle  of  minimum  complementary  energy, 
which  is  a  variational  characterization  of  the  above  problem,  described  by  (2.3)  to 
(2.5),  and  was  first  introduced  by  Hill  ( 1 956),  under  the  assumption  of  strict  convexity 
of  the  nonlinear  potential  U(o,  x).  Thus,  we  define  the  effective  energy  for  the  com¬ 
posite  via  the  relation 

U(a)  =  inf  0(a),  (2.11) 

•eS(rf) 

m 


where 
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0(a)  =  Jfl  U(a,x) dV 

is  the  complementary  energy  functional  of  the  problem  at  hand,  and 

S(d)  =  {<r|<r,77  =  0  in  fl,  and  a.jtij  —  d^nj  on  cfi} 

is  the  set  of  statically  admissible  stresses.  Then,  if  0(6 )  is  assumed  to  be  differentiable, 
it  can  be  readily  shown  that 

eO 

£,7  =  ^=-(<x),  (2.12) 

which  yields  an  effective  constitutive  relation  for  the  composite  in  terms  of  the  average 
variables  a  and  i.  Given  this  connection  between  the  effective  potential  for  the 
composite  0(a)  and  the  effective  stress/strain  relation,  it  makes  sense  to  seek  infor¬ 
mation  on  0(a).  Notice  that,  under  the  above  assumptions,  0(a)  is  convex  (refer  to, 
for  instance.  Appendix  A  of  PCW). 

A  dual  formulation  can  be  given  by  means  of  the  principle  of  minimum  potential 
energy  in  terms  of  the  strain  potential  W(e,x),  which  is  obtained  from  the  stress 
potential  U(a,x)  via  the  Legendre  (Fenchel)  transformation 

lK(s,  x)  =  sup  {a  ■  e—  U(a,x)}. 

Thus,  if  we  define  the  effective  strain  potential  of  the  composite  via 

W(E)  =  inf  IK(e), 

«e*< i) 

where  tV(e)  is  the  pertinent  potential  energy  functional,  and 

K(e)  =  {s|e,7  =  1/2 (u/j  +  ujj)  in  ft,  and  u,  =  e,7.xy  on  dft} 

is  the  set  of  kinematically  admissible  strains  satisfying  a  uniform  strain  boundary 
condition,  we  have  an  effective  stress/strain  relation  for  the  composite,  expressed  by 

div 

6 a  ~  aft*).  (2.15) 

Ofc,7 

where  now  e,  representing  the  average  strain  in  the  composite,  is  equal  to  the  prescribed 
uniform  strain  on  the  boundary,  and  <r,  representing  the  average  stress,  can  be 
“measured”  in  terms  of  the  traction  on  the  boundary.  Notice  that  W(e)  is  also  convex. 

The  above  development  suggests  that  IK(e)  and  0(d)  eaoti  also  be  related  through 
the  Legendre  transformation : 

W(i)  =  sup{<j-E-f?(<f)}.  (2.16) 

■w* 

However,  this  is  certainly  not  true  for  a  general  heterogeneous  material,  since  the 
boundary  conditions  associated  with  the  two  formulations  are  different  (uniform 
traction  versus  uniform  strain),  and  hence  0(d)  and  #'(«)  would  correspond  to  the 
solutions  of  different  problems.  In  fact.  Willis  (1990)  has  shown  that  replacing  the 
uniform  traction  boundary  condition  in  S(d)  by  the  condition  that  the  stresses  have 
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(2.14) 
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mean  value  d  gives  exact  duality.  This  implies  that,  in  general,  the  equality  (  =  ) 
in  (2.16)  must  be  replaced  by  an  inequality  (>).  On  the  other  hand,  it  seems  reasonable 
that  in  the  limit  of  the  microscale  tending  to  zero  (the  homogenized  limit),  the  response 
of  the  composite  would  be  the  same  for  both  types  of  boundary  conditions,  and  thus 
0(d)  and  &(£)  are  expected  to  be  Legendre  duals  in  that  limit.  A  rigorous  proof  of 
this  fact,  however,  would  involve  a  more  rigorous  definition  of  the  homogenized  limit 
than  we  have  utilized,  and  is  beyond  the  scope  of  this  work. 

At  this  point,  a  few  remarks  are  in  order.  First  of  all,  it  should  be  noted  that  the 
three  properties  (2.9),  (2.10)  and  (2.12)  (or  the  corresponding  ones  in  the  dual 
formulation)  hold  true  for  any  heterogeneous  material,  whether  it  is  a  “composite” 
in  the  sense  described  above,  or  not.  However,  0(d)  is  expected  to  represent  the 
effective  properties  of  some  idealized  homogeneous  material,  obtained  as  an  appro¬ 
priately  defined  mathematical  limit  of  a  sequence  of  heterogeneous  materials  with 
vanishingly  small  microscale.  Now,  if  the  microstructure  of  the  composite  is  deter¬ 
ministic  as  in  a  periodic  composite,  0(d)  can  (in  principle)  be  determined  uniquely 
in  terms  of  the  solution  of  a  nonlinear  boundary  problem  on  a  unit  cell  with  periodic 
boundary  conditions  (Marcellini,  1978).  On  the  other  hand,  if  the  microstructure 
of  the  composite  is  random,  usually  only  partial  statistical  information  is  available  in 
the  form  of  the  volume  fractions  of  the  phases,  or,  less  frequently,  some  higher-order 
information  such  as  overall  isotropy  for  the  composite.  It  is  then  not  possible  to 
determine  the  effective  properties  of  any  given  composite  precisely,  and  it  is  essential  to 
reinterpret  0(d)  as  the  set  of  effective  energies  of  a  family  of  composites  with  some 
prescribed  statistics  of  the  microstructure.  In  any  event,  whether  the  composite  is 
periodic  and  its  effective  properties  are  difficult  to  find,  or  random  so  that  its  properties 
are  not  deterministic,  it  makes  sense  to  attempt  to  delimit  the  effective  behavior  of 
composites  by  specifying  bounds  for  0(d)  in  terms  of  some  prescribed  microstructural 
information.  In  some  cases,  as  when  the  bounds  are  too  far  apart  to  be  useful,  it  may 
be  possible  to  identify  a  special  solution,  called  an  estimate,  that  in  some  sense  best 
characterizes  the  properties  of  a  certain  family  of  composites.  In  this  work,  we  will 
only  be  interested  in  the  case  of  composites  for  which  the  volume  fractions  of  the 
constituent  phases  are  specified,  and  that  are,  in  addition,  isotropic  in  an  overall  sense. 
In  the-seqwl.  we  will  attempt  to  specify  bounds  and  estimates  for  the  effective 
properties  of  this  class  of  nonlinear  composites.  ( 


3.  The  New  Variational  Structure 


In  this  section,  we  make  use  of  a  linear  heterogeneous  “comparison”  material,  with 
effective  properties  that  can  be  characterized  in  terms  of  bounds  and  estimates,  to 
obtain  corresponding  bounds  and  estimates  for  the  effective  properties  of  a  nonlinear 
composite.  To  this  end,  we  introduce  the  following  quadratic  potential 


U(c,x)  =  Ix,,,(x)^(<r)  =  — cr;  + 


1 


6/i(x)  '  2k(x) 


(3.1) 


where 
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p(x)  =  £  Z<f>(x)P<,>  >  0.  and  k(x)  =  £  x<',(x)jc,')  >  0, 

f-  I  r-  I 

with  the  /J(,)  and  K{r)  constant,  corresponding  to  a  linear  isotropic  composite  with  the 
same  phase  distribution  (the  same  characteristic  functions)  as  the  nonlinear  composite. 

Then,  if  the  nonlinearity  in  the  potential  of  the  original  composite  is  stronger  than 
quadratic  as  the  norm  of  the  stress  becomes  large,  it  makes  sense  to  define  the  set  of 
functions 

K(r>)  =  sup  {0W(*)-  l/w(*)},  (3.2) 

such  that 

=  t  Z,r>(x)y*'>(fi<'>,K<'>)  =  sup{tf(<7,x)-t/(ff,x)}.  (3.3) 

i  z 

It  follows  that,  for  all  Ktr>  >  0  (r  =  1 , . . . ,  ri)  and  <x,  at  each  xeQ 

U(a,x)  >  U(a,x)-V(M), 

and  hence  that  for  all  fi{'\  K,r)  >  0  (r  =  1 , . . . ,  n),  and  for  every  8 

U(8)>U(8)-V(fi,K),  (3.4) 

where 

0(8)  =  inf  U(o)  (3.5) 

•«SW) 

is  the  effective  potential  of  the  linear  composite,  and 

V(ii ,k)  m  £  cwVw(, 1W,KW), 

r  —  | 


expressed  in  terms  of  the  volume  fractions  of  each  phase, 


J  Cir)  =  J  xw(x)dK 

Thus,  if  we  could  compute  U(8)  for  the  linear  composite  in  terms  of  /i,r)  and  k(r), 
expression  (3.4)  yields  a  family  of  bounds  for  the  effective  potential  of  the  nonlinear 
composite,  U(a),  for  every  choice  of  the  set  of  parameters  /i,r),  jc<,)  >  0.  This  family 
of  bounds  can  be  optimized  by  considering 

0.(8)  =  sup  {( ’j(o)-V(ii,K )}.  (3.6) 


Then,  evidently, 


0(a)  2  0 .(a).  (3.7) 

Usually,  however,  it  is  not  possible  to  find  U(8)  explicitly,  but  instead  bounds  and 
estimates  may  be  available  for  (/(d^If  we  have  a  lower  bound  (such  as  a  Hashin- 
Shtrikman  lower  bound)  U _  (a)  for  U(8),  such  that 


vw 


A'W 


3 


V 


/W A/C 
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U(3)  >  U _  (<r),  (3.8) 

then,  replacing  U ( 3 )  by  U _  (<f  |in  (3^6)  for  t?_  (<r),  still  yields  a  lower  bound  for  0(3) ; 
alternatively,  an  upper  bound  for  0(3)  is  not  useful  in  terms  of  obtaining  an  upper 
bound  for  0(a).  On^the  other  hand,  if  we  only  have  an  estimate  (such  as  a  self- 
consistent  estimate)  0,(3)  for  0(3),  then 

0,(3)  =  sup  {U,(3)-?U 5,  #C)}  (3.9) 

would  provide  only  an  estimate  for  0(3). 

We  note  that  the  prescriptions  (3.6)  and  (3.9)  lead  to  convex  expressions  for  the 
bounds  and  estimates  of  the  effective  potential  of  the  nonlinear  composite,  provided 
that  the  corresponding  bounds  and  estimates  for  the  linear  composite  are  convex, 
which  is  in  turn  guaranteed  assuming  that  p  and  ii  >  0.  This  is  a  desirable  feature, 
because  the  effective  potential  of  the  composite  is  known  to  be  convex,  and  it  follows 
directly  from  the  convexity  of  0(3 )  in  a  (the  supremum  of  a  family  of  convex  functions 
is  convex). 

Although  this  will  not  apply  to  the  present  work,  it  is  possible  that,  in  other  physical 
contexts,  the  nonlinearity  of  the  potential  for  the  composite  will  be  weaker  than 
quadratic.  In  this  case,  upper  bounds  0+(o)  of  the  form  (3.6)  and  estimates  0,(3) 
of  the  form  (3.9)  would  be  obtained  for  0(3),  if  we  replaced  the  suprema  by  infima 
in  the  definitions  (3.2),  and  the  expressions  (3.6)  for  the  bound,  and  (3.9)  for  the 
estimate,  respectively.  In  either  case,  it  is  important  to  note  that  the  given  prescriptions 
for  the  bounds  and  estimates  involve  only  the  evaluation  of  the  extreme  values  of 
some  multidimensional  functions,  assuming  that  the  corresponding  information  is 
available  for  the  linear  comparison  material. 

Before  we  apply  the  general  procedure  developed  in  this  section  to  some  special 
cases,  we  show  that  there  is  no  duality  gap  in  our  procedure.  This  is  unlike  the 
procedure  of  TW,  which  sometimes  leads  to  different  bounds  and  estimates  for  the 
effective  energy  of  the  nonlinear  composite,  depending  on  whether  a  formulation 
based  on  the  principle  of  minimum  complementary  energy,  or  on  the  principle  of 
minimum  potential  energy  is  selected. 

Using  a  procedure  completely  analogous  to  the  above  procedure,  but  starting  with 
the  minimum  potential  energy  formulation,  instead  of  the  minimum  complementary 
energy  formulation,  we  are  led  to  the  following  upper  bound  for  the  effective  potential 

ti'(e) : 

&+(e)  =  i£f>0{^(f)  +  P(A*)},  (3.10) 

where 

fV(e)  =  inf  W(e),  (3.11) 

leXIO 

and 

V(p,K)  =  sup{K'(>,x)-rfr(<,x)}.  (3.12) 

c 

*4 

Notice  that  the  use  of  the  same  notation  for  the  difference  function  V(p,  k)  is  justified. 
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because,  as  the  following  development  shows,  this  present  definition  of  V(p,  k)  is  in 
exact  agreement  with  the  prior  definition  (3.3)  of  V(p,i<).  Thus, 

V(p,k)  =  sup  |  sup  (tr  *  c  —  U  (ff,  x) }  —  W(t,x) 

L  (  Z 

=  sup  jsup{<r-£-  lK(e,x)}  —  U(a,x) 

=  sup  {U{o,x)~  U(a,x)}, 

where  we  have  used  the  fact  that  the  order  of  suprema  may  be  interchanged. 

To  demonstrate  that  there  is  no  duality  gap  in  the  above  procedure,  we  start  with 
the  upper  bound  for  the  effective  strain  potential  of  the  composite 

&(£)  tS  lV+(£).  (3.13) 

Then,  applying  the  Legendre  transformation  to  both  sides  of  this  inequality,  we  get 
(see  Van  Tiel,  1984,  Section  6.3a) 

U(d )  3s  sup  {d‘£- 

where  we  have  made  use  of  (2.16),  and  the  fact  that  both  if'(fi)  and  U  (d)  are  convex. 
Next  note,  from  (3.6),  that 

U.  (d)  =  sup  |  sup  {d-£- W(e)}  —  V(p,  k) 

=  sup  {<*  •  ^^(£‘) +  *)}} 

=  sup  {d •£-&+{£)}, 

where  we  have  made  use  of  some  of  the  properties  of  infima  and  suprema,  and 
assumed  that  duality  in  the  form  of  (2.16)  also  applies  for  the  linear  composite.  Thus, 
we  conclude  that  the  two  bounds  obtained  above  for  the  effective  energy  of  the 
nonlinear  composite  (one  arising  from  the  potential  energy  principle,  and  the  other 
from  the  complementary  energy  principle)  are  equivalent,  and  hence  there  is  no  duality 
gap.  More  generally,  we  can  apply  the  same  ideas  to  the  specialized  versions  of  the 
bounds,  and  again  we  would  obtain  exact  duality  of  the  bounds,  even  for  a  general 
(not  necessarily  a  “composite”)  heterogeneous  material.  For  instance,  if  we  replace 
U(d)  by  its  linear  HS  lower  bound,  and  fF(e)  by  its  linear  HS  upper  bound,  then  the 
exact  duality  of  the  linear  HS  bounds  translates  into  exact  duality  for  the  nonlinear 
bounds.  A  more  formal  development  of  these  ideas,  with  some  additional  results,  is 
given  in  Ponte  Castaneda  (1990). 

4.  Applications 

4. 1 .  The  porous  composite 

In  this  sub-section,  we  apply  the  general  procedure  of  Section  3  to  a  two-phase 
isotropic  composite  with  one  vacuous  phase.  We  take  the  other  phase  to  be  incom¬ 
pressible,  with  potential 
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U°\a)  =  f(at),  (4.1) 

where  /  is  a  strictly  convex  function  of  its  argument  with  stronger  than  quadratic 
growth  as  the  argument  becomes  large.  Then,  the  stress/strain  relation  for  this  phase 
is  given  by 


so  the  the  effective  strain,  e, 
stress  via 


_3  f'(ar) 

z,i  ~  ?  — ~ 
l  a. 


Sii> 


(4.2) 


'Oe  •  e  (e  is  the  deviator  of  a),  is  related  to  the  effective 


«,  =  if  (<*,)• 

On  the  other  hand,  the  potential  of  the  vacuous  phase  is  given  by 


Injsw\ 

where 


Now,  if  we  let 


<5oM  = 


U[2)(a)  =  S0(o.),+  &o(<T~) 

0  if  .X  =  0 
oo  otherwise. 


(4.3) 


and 

we  have  that 


(7(2)(ff)  =  <50(ff,^ 

ss  \af'(a)—f(a). 


(4.4) 


(4.5) 


(4.6) 


where  <r  is  some  function  of  p(l)  determined  by  the  optimization  problem  (3.2),  and 
also  that 


K,2,(/j,2))  =  0.  (4.7) 

It  follows  from  (3.6)  that 

0.(a)  =  sup  {£(<7)  — cll)r"Vl))},  (4.8) 

(i’">  o 

provides  a  general  lower  bound  for  0(a)  given  0(a).  As  discussed  previously,  in 
general,  we  do  not  know  U(a)  precisely.  Here,  we  will  make  use  of  the  Voigt  and 
lower  Hashin-Shtrikman  (HS)  bounds  for  0(a)  to  obtain  corresponding  bounds  for 
U(a).  Additionally,  we  will  provide^  self-consistent  (SC)  estimate  for  0(a)  in  terms 
of  the  well-known  SC  estimate  for  0(a). 


4.1.1.  Voigt  bound.  For  the  linear  composite  with  one  vacuous  phase,  it  is  known 
that 
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U(6)>^-dl 

where  {iv  —  c(n/i(l).  Therefore,  from  relation  (4.8),  we  have  that 

U{a)^c['}  sup  {7-TrT^-Kt,>(M,1,)j 
^'">o  I6/1  3 

-c“m 

where 

s  =  djc"\ 
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(4.9) 

(4.10) 


and  where  we  have  made  use  of  the  result  of  the  Appendix,  and  hence  we  must  further 
assume  that  F(x)  =  f(s)  is  a  convex  function  of  x  =  s:  >  0.  The  result  expressed  by 
(4.9)  and  (4.10)  is  precisely  the  nonlinear  Voigt  bound,  which  could  alternatively  be 
obtained  from  the  principle  of  minimum  potential  energy  by  assuming  a  uniform 
strain  field  throughout  the  composite,  and  dualizing  the  result.  Thus,  at  least  in  this 
simple  case,  the  new  procedure  reproduces  the  “right”  result,  exactly. 


4.1.2.  Hashin-Shtrikman  bound.  For  linear  isotropic  composites,  Hashin  and 
Shtrikman  (1962)  found  upper  and  lower  bounds  that  are  tighter  than  the  Voigt/ 
Reuss  bounds.  The  lower  bound  for  our  particular  example  specializes  to 


U  (S)  > 


fy'lHS 


d;  + 


where 


fi  hs  — 


c"y,) 


and  — 


4  c1 


,<n 


1  +  5C<:>  HS  ~  3  pcc- 


Tf1 


(I) 


Then,  a  bound  of  the  form  (4.9)  applies,  where  now 


J  =  -7»y(l  +  5Ci;i)<T,:  +  4Ci;:!<?™, 


(4.11) 


and  again  we  have  made  use  of  the  result  of  the  Appendix.  Because  this  bound  was 
derived  from  the  linear  HS  bound,  we  refer  to  it  as  the  nonlinear  HS  bound.  Note 
that,  unlike  the  Voigt  bound,  it  predicts  overall  compressibility  for  the  composite,  as 
expected  physically. 


4.1.3.  Self-consistent  estimate.  For  linear  isotropic  composites.  Budiansky  (1965) 
and  Hill  (1965)  provided  the  so-cailed  SC  estimates  for  the  effective  moduli.  In  this 
particular  case,  these  estimates  specialize  to 


U{6) 


6  fisc 


dt  + 


2k, 


Gfn  ♦ 


where 


\2 


i  ow'Oe'^ 


a„Jt  f''v"•',  »*V,CUS^ 


P. 


/ 
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fisc  — 


(1  — 2f<2))  „„  .  4f'"  (1  — 2c!2>) 

(l-c,:,/3)P  ’  3nd  Ksc  3  c‘2>  (1  —  c(2)/3) ^  ‘ 


Then 

0(a)  »  c(I>/(j), 

where  now 


s  = 


,to 


1  —  r‘:,/3 
1  —  2ct:i 


9  c,:i 

ff<+ 4  7ITa- 


(4.12) 


and  once  again  we  have  made  use  of  the  Appendix.  We  refer  to  this  result  as  the 
nonlinear  SC  estimate. 


4.2.  The  two-phase  incompressible  composite 

In  this  sub-section,  we  consider  a  general  two-phase  composite  with  isotropic, 
incompressible  and  nonlinear  phases  such  that 

Uv'(a)  =r\ot),  (r  =  1,2),  (4.13) 

where  the  /'"  satisfy  the  same  convexity  conditions  of  the  previous  section.  Then  if 
we  let 


(4.14) 


we  have  that 

(4.15) 

where  the  cr1'’  are  some  functions  of  the  fi{r)  that  are  determined  by  the  solution  of  the 
optimization  problem  (3.2).  It  follows  that 


is  a  general  lower  bound  for  0(a),  given  U(a).  As  in  the  previous  sub-section,  we  use 
this  result  to  obtain  the  Voigt  bound,  a  HS  bound  and  a  SC  estimate  for  the  effective 
potential  of  the  nonlinear  composite  in  terms  of  the  corresponding  results  for  the 
linear  composite. 


4.2.1.  Voigt  bound.  In  this  case,  it  is  more  convenient  to  consider  the  dual  formu¬ 
lation.  For  the  iinear  composite,  we  have  that 

\V(i)  sz  ]pye;+60(£m),  (4.17) 

where  py  =  c'"p"'  +  c,1>p'Z).  Then,  it  can  easily  be  shown  hat 

W(£)  s:  c",M'<"(e)  +  c(J'W'a,(«).  (418) 

where  w  e  have  made  use  of  the  appropriate  specialization  of  (3.10),  and  a  dual  version 
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of  the  result  of  the  Appendix.  As  in  the  previous  example,  this  result  agrees  exactly 
with  the  nonlinear  Voigt  bound. 

4.2.2.  Hashin-Shtrikman  bounds.  For  the  linear  composite,  we  have 

U(a)  (4.19) 

o  l*HS 

where 

c‘  V '  W 31  +  9/i )  +  c‘  V 3)  (6/i‘ 1 '  +  9/i ) 
hHS  ~  c,‘,(6/i'3’  +  9ji)  +  c'3W,’+9/i) 

and  p  ■  x  {/i(l\/i(2)}  =  piU  (by  assumption).  Then,  it  follows  that 

#(<?)>  sup  ef,--c<l>K",(^<,>)-c‘=)  K'2*(m(2>)  }  (4.20) 

o  (6 pHS  J 

but,  unfortunately,  no  further  simplification  is  possible  in  general,  due  to  the  coupling 
of  pil)  and  /i131  in  the  term  involving  pHS.  Later  we  will  apply  this  result  to  a  special 
case  where  /,l,,/tJ)  have  a  simple  power-law  form. 

4.2.3.  Self-consistent  estimates.  The  results  for  the  SC  estimates  have  the  same 
forms  (4.19)  and  (4.20)  as  the  HS  bounds,  but  with  pHS  replaced  by  pSc*  given  by  the 
positive  root  of  the  expression 

3/iJc  +  {(2  -  5c‘  ‘V  "  +  (2  -  5c'3  V:,}/<sc  —  V:i  =  0.  (4.21) 

As  was  the  case  for  the  HS  bound,  these  results  cannot  be  simplified  further  without 
specifying  in  more  detail  the  constitutive  behavior  of  the  phase  materials. 


4.3-  The  composite  reinforced  by  rigid  inclusions 

This  is  a  special  case  of  the  previous  material  with  P'\o,)  =  0,  corresponding  to 
the  case  where  phase  #1  is  rigid.  Then,  the  choice  Un)(a)  =  0  leads  to  ^"’(/i"1)  =  0. 
and  we  have  the  following  general  lower  bound  for  the  effective  energy  of  the  com¬ 
posite  : 

C(a )^  sup  {(/(it)  —  c':i  F<3l(/ii31)}-  (4.22) 

li'  *'  >  0 

In  this  case,  however,  it  is  clear  that  the  lower  bounds  are  trivial,  and  we  will  only 
be  able  to  give  estimates  for  the  effective  energy  of  the  composite.  We  will  consider 
the  standard  SC  estimate,  as  well  as  a  differential  self-consistent  (DSC)  estimate. 

4.3.1.  Self-consistent  estimate.  For  the  linear  composite,  it  is  known  that 

U(d )  ^  -t—  a;,  (4.23) 

opsc 


where 
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,',(21 


/2jc  — 


(i-k"') 


&"<  D- 


Using  this  result  in  the  general  expression  (4.22)  induces  a  nonlinear  SC  estimate  for 
the  effective  energy 


where 


0(3)  a  cr-'f2>(s). 


(4.24) 


(4.25) 


and  once  again  we  have  made  use  of  the  result  in  the  Appendix. 


4.3.2.  Differential  self-consistent  estimate.  The  DSC  estimates  of  Boucher  (1974), 
and  McLaughlin  (197,4^  for  linear  elastic  composites  specialize  in  this  case  to 
expression  (4.23),  with  psc  replaced  by 

p'J> 

Vdsc  -  ^TTiyn  •  ~  (4.26) 

This  induces  a  nonlinear  DSC  estimate  via  (4.22)  that  reduces  to  an  expression  similar 
*o  (4.24),  where  now 

s  =  (cl2>)34df.  (4.27) 


5.  Results  for  Power-law  Materials 

In  this  section,  we  specialize  further  the  calculations  of  the  previous  section  by 
taking  the  constitutive  behavior  of  the  phases  to  be  governed  by  a  power  law  relation 


This  class  of  functions  clearly  satisfies  all  the  assumptions  invoked  in  Sections  3  and 
4,  including  the  convexity  assumption  of  the  Appendix.  Additionally,  we  compare  the 
new  results  for  the  bounds  and  estimates  with  previously  available  results. 

5. 1 .  The  porous  material 

In  this  case,  our  results  for  the  Voigt  bound,  the  HS  bound  and  the  SC  estimate  all 
take  the  simple  form  of  expression  (4.9),  with  s  given  by  (4.10),  (4.11)  and  (4.12), 
respectively.  We  compare  these  results  with  the  results  of  PCW  for  the  same  material. 
We  note,  however,  that  Willis  (1989)  has  given  a  HS  bound  and  several  SC  estimates 
for  the  general  case  considered  in  Section  4.1,  but  the  form  of  his  results  is  more 
complicated  than  our  new  results.  For  this  reason,  and  because  we  expect  the  com¬ 
parison  of  the  power-law  results  to  be  fairly  representative,  we  do  not  make  a  more 
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general  comparison  between  our  new  results  and  those  of  Willis  (1989).  There 
are  two  domains  of  particular  interest:  the  low-triaxiality  range  corresponding  to 
a)  =  l<fm/<i,l  «  1 ,  and  the  high-triaxiality  range  corresponding  to  id  »  1 . 

For  low-triaxiality,  all  the  results  for  the  bounds  and  estimates  take  the  form 

<5-2) 


and  it  suffices  to  compare  all  the  corresponding  values  of  the  low-triaxiality  modulus 
ji ,  as  functions  of  pn),  ci2)  and  n.  We  give  the  results  below,  where  we  identify  the 
source  of  the  results  in  parentheses : 


The  last  expression  corresponds  to  a  dilute  concentration  of  voids,  and  was  carried 
out  by  DH,  who  made  use  of  the  results  of  Budiansky  el  al.  (1982).  The  values  of 
/*,  k  and  x*  as  functions  of  n  are  taken  from  that  reference. 

In  the  linear  limit  (n  =  1),  all  the  above  results  agree  with  the  well-known  results 
from  the  linear  theory.  For  small  volume  fractions  of  the  vacuous  phase  (c(21).  the 
HS,  SC  and  dilute  results  agree  to  first  order,  and  for  larger  volume  fractions  the 
dilute  result  lies  between  the  HS  upper  bound  and  the  SC  estimate.  Figure  1  gives  the 
results  for  the  bounds  and  estimates  of  the  effective  low-triaxiality  moduli  fi  as 
functions  of  cu\  for  two  distinct  nonlinear  cases  (n  =  3  and  10).  These  results  are 
similar  to  the  linear  results,  with  the  HS  upper  bounds  lying  below  the  Voigt  bound, 
and  the  SC  estimates  close  to  the  dilute  results  for  moderate  volume  fractions,  but 
approaching  the  same  percolation  limit  (c(2)  =  1/2)  as  the  corresponding  linear  results. 
Comparing  the  new  bound  and  estimate  with  those  of  PCW,  we  observe  that  the  new 
results  lie  slightly  below  the  old  results.  For  the  bound,  this  has  the  implication  that 
the  new  bound  is  an  improvement  on  the  old  bound,  at  least  in  the  low-triaxiality 
range. 
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',Br0U"ds  an?  “"males  for  the  low-tnaxialitv  modulus  of  Ihe  porous  material  as  functions  of  th 
olume  fraction  of  voids.  The  short-dash  line  corresponds  to  the  Voigt  bound;  the  continuous  line 

rW  Tu  and  SC  Ktimale ;  the  lon?-dastl  line* correspond  to  the  HS  uppe 
and  SC  estimate  of  PCW  .and  the  long/short-dash  line  corresponds  to  the  dilute  calculation  of  DH 
Case  (a)  is  for  n  =  3,  and  (b)  for  n  =  10. 
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For  high  triaxiality,  all  the  results  for  the  bounds  and  estimates  take  the  form 
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/  (27)'* 

— 

2  m 
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(5.4) 


and  thus  we  express  all  the  results  in  terms  of  the  high-triaxiality  moduli  y,  so  that 


/ 1' 

7  HS 


I, ID 


+  I)  2/i  1 


In  the  linear  limit  (n  =  I),  the  expressions  for  the  HS  bounds  and  SC  estimates  reduce 
to  the  linear  results,  and,  additionally,  they  agree  in  the  dilute  limit  (y  ~  (ci:,)_1). 
Figure  2  depicts  plots  of  the  HS  bounds  and  SC  and  dilute  estimates  of  the  high- 
triaxiality  modulus  appropriately  normalized  ((c,:>)("+ 1,1  *”f)  versus  the  volume  frac¬ 
tion  of  the  void  phase  (ci:>)  for  two  values  of  the  nonlinearity  parameter  (n  =  3  and 
10).  By  comparison  with  the  results  of  PCW,  we  find  that  the  new  results  provide  a 
significant  improvement  over  the  old  results,  since  the  new  bounds  lie  significantly 
below  the  old  bounds,  and  are  hence  tighter.  Similarly,  the  old  SC  estimate  violates 
the  new  bound,  and  must  be  discarded  in  favor  of  the  new  SC  estimate.  The  dilute 
result  of  DH  lies  below  the  HS  bound  for  small  volume  fractions  of  the  void  phase 
(Vo  ~  (c(2,)HI  whereas  yHS  ~  (c(:,)_("+ 1,;>),  but  it  is  clear  from  the  plots  that  the 
range  of  validity  of  the  dilute  result  is  severely  limited  for  larger  values  of  n.  On  the 
other  hand,  we  should  emphasize  that  the  new  expressions  for  y  are  not  valid  for  small 
values  of  c<:).  This  is  because  the  original  expression  for  U(a),  from  which  they  derive, 
is  indeterminate  in  the  limit  as  co  -*  oo  and  c<:’  -*  0.  A  better  comparison  of  the  new 
results  with  the  dilute  results  is  accomplished  by  making  use  of  the  original  expression 
for  0(a)  in  the  form 


Fig.  2.  Bounds  and  estimates  for  the  high-triaxiality  modulus  of  the  porous  material  as  functions  of  the 
volume  fraction  of  voids.  The  continuous  lines  correspond  to  the  new  HS  upper  bound  and  SC  estimate; 
the  long-dash  lines  correspond  to  the  HS  upper  bound  and  SC  estimate  of  PCW ;  and  the  short-dash  line 
corresponds  to  the  dilute  calculation  of  DH.  Case  (a)  is  for  n  =  3.  and  (b)  for  n  =  10. 
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(2/i*  r 


<r/ 


where  now  ft*  is  a  function  of  c(2),  n  and  to.  Thus,  we  have  that 


[l+c(J,(5+jar)],',+  ')'2"’ 


(5.6) 


(5.7) 


and 


/<"’  [1  +  c,: 21  (n+l)/(co,  n)]1'*’  K  * 

where  f(co,n)  is  taken  from  the  work  DH. 

Written  in  this  form,  it  is  clear  that  both  expressions  for  the  effective  energy  of  the 
composite  are  indeterminate  in  the  limit  as  co  -*•  co  and  r!2>  -» 0.  The  first  expression, 
however,  is  general,  whereas  the  second  assumes  that  c(2)  «  1,  and  its  range  of  validity 
is  not  known.  Figure  3  shows  a  comparison  of  these  two  results  as  functions  of  c o  for 
two  values  of  n  (3  and  10)  and  three  values  of  c(2)  (0.1,  0.01  and  0.001).  It  is  apparent 
that  the  dilute  approximation  of  DH  is  acceptable  (although  this  is  not  a  proof  that 
it  is  correct)  for  very  small  values  of  c,:>  in  the  sense  that  it  does  not  violate  the  new 
bound,  but  unacceptable  for  values  of  c(2)  in  the  order  of  1  to  10%,  or  larger,  depending 
on  the  specific  value  of  n.  On  the  other  hand,  it  appears  that  the  new  bound  could 
conceivably  be  subject  to  improvement  for  very  small  values  of  c<:>. 


5.2.  The  two-phase  incompressible  composite 

In  the  case  when /l0  and /(2)  have  the  same  form  of  (5.1),  but  with  different  moduli 
piU  and  /i<2),  respectively,  all  the  results  of  Section  4.2  take  the  form 


U(S)  = 


1 


3  \n  +  1/  (2/3) 


1 


(5.9) 


and  it  suffices  to  compare  the  bounds  and  estimates  for  the  effective  modulus  ft.  The 
Voigt  bound  is  obtained  from  (4.18),  and  is  given  by 


Hr 

„ii> 


=  c(l)  +c' 


(2)  ! 


,(2. 


(5.10) 


The  new  HS  bound  is  obtained  by  solving  the  optimization  problem  given  in  (4.20) ; 
the  result  is 


Hhs 


(5.11) 


Mf  tr 


fcffrv  SCINCH 


where  S(  n  and  S,2)  satisfy  the  relations 
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Fig.  3.  Bounds  and  dilute  estimates  for  the  high-triaxiality  modulus  of  the  porous  material  as  functions 
of  the  triaxiality  for  three  different  volume  fractions  of  voids.  The  continuous  line  corresponds  to  the  new 
HS  upper  bound,  and  the  dash  line  corresponds  to  the  dilute  calculation  of  DH.  Case  (a)  is  for  n  =  3.  and 

(b)  for/t  =  10. 
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[5-(2  +  3c«2>)5(:)]3?nir= 


and 


15c«"S(Ju  =  6C(I>5(:,+  -jfy  (3  +  c<2)(2  — 5S(2,)][5  -(2  +  3c,2))5,2)],  (5.12) 


subject  to  the  restriction  that 


The  SC  estimate  usc  is  obtained  from  (4.20),  and  the  result  has  the  same  form  as 
(5.11),  but  now  S(i/lnd  5, 2)  satisfy 


c'^fM  +  C1 


(2)  ^ 


(I) 


_  C  2  _ 

-12)  15  (2)  ~ 


p}  " 

fisc 


and 


3^+(2-5  c(") 


=  co 

J  fi{" 


(2) 


6^+(2-5c-")  +  (2-5c‘2>)^ 


L 


P 'Jfi 


fisc 

I>  ^l2)» 


(5.13) 


where 


,-.l  2) 


1 


and  fisc  is  the  positive  root  of  (4.21).  Additionally,  for  this  problem  we  have  a  non¬ 
trivial  Reuss  lower  bound  (if  /t<:>  >  0)  given  by 


CU»+C(2) 


—  (In) 


(5.14) 


In  the  linear  limit  (rt  =  1),  the  above  results  reduce  to  the  well-known  results  from 
the  linear  theory  and  additionally  there  is  a  HS  lower  bound  for  fi.  The  HS  bounds 
are  nested  within  the  Reuss/Voigt  bounds,  and  the  SC  estimate  lies  within  the  HS 
bounds,  agreeing  with  the  HS  upper  bound  for  small  volume  fractions  of  phase 
#2  (c<2)).  and  with  the  HS  lower  bound  for  small  volume  fractions  of  phase  #1 
(c'1’  =  1  —  c<2))-  As  the  ratio  becomes  smaller,  the  bounds  spread  out  until 

eventually  (when  il2)/n(l)  =  0)  the  lower  bounds  become  trivial  and  the  SC  estimate 
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reaches  a  percolation  limit  at  c(2)  =  3/5.  Figure  4  depicts  the  corresponding  results 
(except  the  HS  lower  bound)  for  the  two  usual  values  of  n  (3  and  10).  By  comparison 
with  the  linear  results  (not  shown),  they  roughly  show  the  same  trends.  However,  the 
following  comparisons  can  be  established:  the  Voigt  bound  is  unaffected;  the  SC 
estimate  is  shifted  up  for  low  values  of  cl2\  and  down  slightly  for  low  values  of  cm, 
with  larger  values  of  n;  and  the  Reuss  bound  is  shifted  down  significantly  with 
larger  values  of  n  to  the  point  of  being  nearly  without  practical  value  for  large  n. 
Extrapolating  from  the  linear  theory,  and  according  to  the  new  results,  we  expect  the 
SC  estimates  to  be  good  at  predicting  the  effective  moduli  of  the  composite,  if 
is  not  too  small. 


5.3.  The  composite  reinforced  by  rigid  inclusions 

In  this  case,  as  in  the  previous  one,  all  the  results  take  the  form  of  relation  (5.1), 
where  now  we  no  longer  have  bounds  for  fi,  but  only  estimates,  including  a  SC  as 
well  as  a  DSC  estimate. 

The  SC  estimate  is  obtained  from  expressions  (4.24)  and  (4.25)  and  is  given  by 


fisc  (c'2')"-"2' 


(5.15) 


The  corresponding  result  from  PCW  is  given  by 


gs-d-^4  ■»)*•  (5-,6) 

Our  DSC  estimate  is  obtained  from  the  corresponding  linear  estimate  via  expression 
(4.24),  together  with  (4.27)  and  is  given  by 


fipsc  1 

^(2)  —  (c<2)yj<!+7);4«  • 


(5.17) 


Duva  (1984)  has  also  provided  a  DSC  estimate,  based  on  computations  for  an 
isolated  rigid  inclusion  in  an  infinite  power-law  viscous  matrix.  The  result  of  his 
calculation  can  also  be  represented  in  the  form  of  (5.1),  with 


fipsc _ 1 

—  ^c(2)^?l*l  <i  > 


(5.18) 


where  g(n)  is  given  in  the  above  reference.  From  the  plot  of  this  function,  we  obtain 
that  0(1)  =  5/2, 0(3)  ft  3.3  and  0(10)  ft  6.6. 

As  before  our  new  results  reduce  to  the  linear  results  in  the  limit  as  n  approaches 
unity.  Figure  5  gives  results  for  the  above  estimates  of  the  effective  modulus  for  the 
two  usual  values  of  n  (3  and  10).  Note  that  the  plot  is  given  in  terms  of  as  a 
function  of  ci:\  so  that  p,2)/fi  tends  to  zero  as  c<2)  tends  to  unity.  Our  SC  estimate  is 
slightly  lower  than  the  prior  result  of  PCW ;  this  is  probably  due  to  the  fact  that  if  we 
had  an  upper  bound,  it  would  also  appear  lower  in  our  plot.  On  the  other  hand,  our 
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Fig.  4.  Bounds  and  estimates  Tor  the  effective  modulus  of  the  two-phase  incompressible  composite  as 
functions  of  the  volume  fractions  of  phase  #  2  for  three  different  ratios  of  the  moduli  of  the  two  phases. 
The  continuous  lines  correspond  to  the  new  HS  upper  bound  and  SC  estimate,  and  the  dotted  lines 
correspond  to  the  Voigt/Reuss  bounds.  Case  (a)  is  for  n  =  3,  and  (b)  for  n  =  10. 


Fig.  5.  Estimates  of  the  effective  modulus  of  the  rigidly  reinforced  composite  as  functions  of  the  volume 
fraction  of  the  rigid  phase.  The  continuous  lines  correspond  to  the  new  SC  and  DSC  estimates,  and  the 
dash  lines  correspond  to  the  SC  estimate  of  PCW,  and  the  DSC  estimate  of  Duva  (1984).  Case  (a)  is  for 

n  =  3.  and  (b)  for  n  =  10. 
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DSC  estimate  lies  somewhat  above  the  DSC  of  Duva  (1984).  This  is  not  too  surprising 
since  we  expect  our  procedure  in  general  to  fin  ilimuluthc  energy  of  the  composite 
&(£).  For  the  linear  case,  it  has  been  shown  (Milton,  1985;  Avellaneda,  1987) 
that  the  DSC  estimates  can  be  attained  by  particular  microstructures.  Then,  we  could 
expect  our  procedure  to  give  an  upper  bound  for  the  actual  #(£)  of  the  corresponding 
nonlinear  material.  Thus,  the  present  example  could  be  used  to  assess  the  efficiency 
of  our  method  irjestimating  the  exact  effective  properties  of  nonlinear  composites  with 
deterministic  rrucroztructure,  such  as  periodic  composites.  This  idea  will  be  tested 
elsewhere,  but  if  the  indications  of  the  present  example  are  upheld,  we  would  conclude 
that  our  prescription  performs  rather  well  given  its  simplicity. 

6.  Concluding  Rem  vrks 

The  main  contribution  of  this  work  is  the  establishment  of  a  new  variational 
structure  that  allows  the  estimation  of  the  effective  properties  of  nonlinear  composites 
in  terms  of  the  corresponding  properties  for  linear  composites  with  the  same  micro- 
structural  distribution  of  phases.  Assuming  that  the  exact  effective  energy  density  is 
available  for  the  linear  composite,  and  depending  on  the  growth  conditions  of  the 
phase  potentials,  the  new  estimate  will  either  be  an  upper  bound,  or  a  lower  bound 
for  the  actual  effective  energy  density  of  the  nonlinear  composite.  Alternatively,  if 
only  an  estimate  or  a  bound  of  the  right  type  (an  upper  bound  if  the  estimate  is  an 
upper  bound,  or  vice  versa)  is  available  for  the  linear  composite,  the  new  estimate  for 
the  nonlinear  composite  will  also  be  either  only  an  estimate,  or  a  bound  (of  the  same 
type).  In  this  respect,  the  present  variational  structure  has  the  same  limitation  as  the 
structure  proposed  by  TW  in  that  only  one-sided  bounds  will  result  in  general. 

In  the  context  of  specific  results,  we  note  that  the  nonlinear  bounds,  corresponding 
to  the  linear  Voigt  bounds,  proposed  by  the  new  prescription  turn  out  to  be  precisely 
the  nonlinear  Voigt  bounds  obtained  directly  from  the  principle  of  minimum  potential 
energy.  On  the  other  hand,  the  nonlinear  bounds  corresponding  to  the  Hashin- 
Shtrikman  bounds  for  the  linear  isotropic  composite  turn  out  to  be  superior  in  some 
cases  to  the  Talbot-Willis  bounds  for  the  same  material,  and  identical  in  other  cases 
(see  Ponte  Castaneda,  1990,  for  an  example).  Finally,  whereas  the  Talbot-Willis 
structure  leads  to  some  ambiguity  in  the  prescription  of  self-consistent  estimates,  the 
new  structure  leads  to  a  unique  prescription  for  a  self-consistent  estimate,  as  well  as 
to  a  straightforward  generalization  of  other  types  of  estimates,  including  differential 
and  dilute  estimates.  The  main  advantages  of  the  new  structure  are  the  simplicity  of 
its  implementation  and  the  generality  of  its  potential  applications.  Clearly,  these  are 
features  that  could  make  the  proposed  structure  of  great  practical,  as  well  as  theoretical 
value.  •A-.  iiuui.  uliunwt  derivation  of  this  structure  is  given  in  Ponte  Castaneda 
(1990),  and  some  of  the  potential  applications  to  other  types  of  composites  will  be 
addressed  elsewhere. 
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Exact  Duality  Between  /  and  V"’ 
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Here  we  demonstrate  that 


sup 

^"’>0 


given  that 
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K‘l>(p,,,)  =  sup|^iTJ:-/(i)|.  (A2) 


proof :  Let  x  =  r,  and  assume  that  /"(*)  **  /(s)  is  a  convex  function  of  its  argument.  Then,  if 
we  define  G(y)  via  the  Legendre  transform 

G(y)  =  sup  {.xy-.F(.x)},  (A3) 

x  >  0 

we  have  by  Fenchel  duality  (Van  Tiel,  1984,  Corollary  to  Section  6.1  la)  that 

F(.x)  =  sup  {.xy-G(y)}.  (A4) 

y  >  o 

Now,  if  we  let  y  =  1/6/i1",  (A2)  and  (A3)  imply  that 

G(y)  -  F"V>),  (A5) 


and  therefore  it  follows  from  (A4)  that  (1)  holds. 

>  This  result  is  a  special  case  of  a  more  general  result  derived  in  Ponte  Castaneda  (1990). 
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Abstract 

A  new  variational  procedure  recently  developed  by  Ponte  Castaneda  [1]  for  estimating  and 
bounding  the  effective  properties  of  nonlinear  composite  materials  is  applied  to  an  incompressible 
isotropic  matrix  containing  an  isotropic  distribution  of  voids.  The  uniaxial  stress-strain  behavior  of 
the  matrix  is  linear  up  to  the  yield  stress,  and  "linear  plus  power-hardening"  for  stresses  exceeding 
the  yield  level.  The  solution  procedure  exploits  a  useful  "comparison"  between  the  effective 
potential  of  the  nonlinear  porous  material  and  that  of  an  appropriately  chosen  inhomogeneous, 
linear  material.  It  is  important  to  note  that  the  proposed  procedure  is  different  from  the  procedure 
advocated  by  Talbot  &  Willis  [2],  which  makes  use  of  a  linear,  but  homogeneous  comparison 
material  to  accomplish  the  same  goal.  Results  are  given  in  the  form  of  a  rigorous  lower  bound  and 
a  self-consistent  estimate  for  the  effective  potential  of  the  porous  material.  The  results  are  compared 
with  the  results  of  Willis  [3]  for  the  same  problem. 

Introduction 

The  theory  of  linear  composites  is  fairly  well  developed,  including  different  approaches  to 
the  problem.  Thus,  exact  estimates  have  been  determined  for  the  effective  properties  of  ad  hoc 
models  for  composites;  rigorous  variational  bounds  have  been  given  for  the  properties  of  random 
composites;  and  explicit  formulae  have  been  given  for  the  properties  of  periodic  composites. 
Appropriate,  but  by  no  means  exhaustive,  references  dealing  with  the  linear  theories  are  provided 
by  [4-6].  By  contrast,  the  theory  of  nonlinear  composites  is  not  very  well  developed,  and  most  of 
the  results  are  based  on  ad  hoc  models,  such  as  dilute  and  self-consistent  models.  For  instance, 
dilute  estimates  are  given  in  [7]  for  the  effective  properties  of  a  nonlinear  porous  material.  The  first 
contribution  dealing  with  the  calculation  of  rigorous  bounds  for  the  effective  properties  of 
nonlinear  composites  is  provided  by  the  work  of  Talbot  and  Willis  [2].  These  authors  extended  the 
well-known  Hashin-Shtrikman  [8]  variational  principles  to  include  nonlinear  constitutive  behavior, 
and  their  methods  have  been  applied  to  a  number  of  examples  in  different  physical  contexts.  For 
example,  bounds  and  estimates  for  the  effective  properties  of  nonlinearly  viscous  (or  deformation- 
theory  plastic)  materials  are  determined  in  [3]  and  [9].  A  new  variational  procedure  for  estimating 
the  effective  properties  of  nonlinear  composite  materials  was  developed  recently  by  Ponte 
Castaneda  [1].  This  procedure  can  be  used  directly  to  obtain  bounds  and  estimates  for  the  effective 
properties  of  nonlinear  composites  from  any  corresponding  bounds  and  estimates  that  may  be 
available  for  the  effective  properties  of  linear  composites  with  the  same  distribution  of  phases.  The 
main  advantages  of  the  new  procedure  over  the  procedure  given  by  [2]  are  the  simplicity  of  its 
implementation,  the  generality  of  its  potential  applications  and  the  strength  of  the  results.  In 
reference  [1],  the  general  procedure  is  also  applied  to  an  incompressible  isotropic  matrix  containing 
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an  isotropic  distribution  of  voids,  and  specific  results  in  the  form  of  a  Hashin-Shtrikman  bound 
and  a  self-consistent  estimate  are  provided  for  pure-power  law  behavior  for  the  matrix  material. 
The  results  for  the  bounds  are  found  to  be  stronger  than  the  corresponding  bounds  given  in  [9].  In 
this  work,  we  consider  the  application  of  the  general  result  for  a  porous  material  for  the  special 
case  of  stress-strain  behavior  for  the  matrix  which  is  linear  up  to  the  yield  stress  and  linear  plus 
power-hardening  for  stresses  exceeding  the  yield  level.  The  new  results  are  discussed  and 
compared  with  the  corresponding  results  of  Willis  [3]. 


New  Nonlinear  Hashin-Shtrikman  Bounds  and  Self-Consistent  Estimates 

We  are  interested  in  estimating  the  effective  properties  of  composites  with  nonlinear 
material  behavior.  By  a  "composite"  we  mean  an  idealized  material  that  corresponds  to  the  limit  of 
a  sequence  of  heterogeneous  materials  with  two  distinct  length  scales:  one  microscopic  / 
corresponding  to  the  size  of  the  heterogeneity,  and  one  macroscopic  L  corresponding  to  the  size  of 
the  specimen  of  interest  and  the  scale  of  variation  of  the  boundary  conditions.  The  effective 
behavior  of  the  composite  is  then  obtained  by  considering  the  limiting  behavior  of  the  sequence  of 
materials  as  the  ratio  of  scales  e  =  UL  tends  to  zero. 

Consider  an  n-phase  composite  occupying  a  domain  Q  (normalized  to  have  unit  volume), 
with  each  phase  occupying  a  sub-domain  Q(r)  (r  =  1,  2,...,  n ),  and  let  the  stress  potential, 

f/(a,x),  be  expressed  in  terms  of  the  n  homogeneous  phase  potentials,  £/(r)(a),  via 


where 


U(o,x)  =  '£X<r\x)U<rXo), 

r=  1 


(1) 


X<r\x)  = 


(2) 


if  x  €  Q(r> 
otherwise 

is  the  indicator  function  of  phase  r.  The  phases  are  assumed  to  be  isotropic,  so  that  the  potentials 


£/(,)(< 7)  depend  on  the  stress  O  only  through  its  three  principal  invariants.  Here,  we  will  further 
assume  that  the  dependence  is  only  through  two  of  the  invariants,  namely,  the  mean  stress. 


1  h 

om  =  —  tra,  and  the  effective  stress,  <7  = ,  —  S  •  S ,  where  S  is  the  deviator  of  a. 

3  V  2 

The  effective  behavior  of  the  composite  material  is  defined,  following  [10],  in  terms  of  the 
effective  energy,  (7(a),  that  arises  due  to  the  uniform  traction  boundary  condition 

atjni  =  G-  n.j,  x  e  dQ,  (3) 

where  dQ  denotes  the  boundary  of  the  composite,  n  is  its  unit  outward  normal,  and  ct  is  a  given 

constant  symmetric  tensor.  The  average  stress  in  the  composite  is  then  precisely  a  and  it  follows 
from  the  principle  of  minimum  complementary  energy  that 


U(a)  =  min  U(a), 

oeS(a) 


(4) 


where 


C/(o)=f  U(o,x)dV 
Jo 
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is  the  complementary  energy  functional  of  the  problem,  and 

5(a)  =  {al <Tyj  =  0  in  £2,  and  a^n ;  =  aiyny  on  dQ) 

is  the  set  of  statically  admissible  stresses.  Thus,  if  e  denotes  the  average  strain  over  the  composite, 
it  can  be  readily  shown  that 


-  dU  ,-V 
V 


(5) 


which  yields  an  effective  constitutive  relation  for  the  composite  in  terms  of  the  average  variables  a 
and  e.  Given  this  connection  between  the  effective  potential  for  the  composite  (7(a)  and  the 


effective  stress/strain  relation,  it  makes  sense  to  seek  information  on  (7(a).  Notice  that  (/(a)  is 
convex. 

Next,  following  [1],  we  make  use  of  a  linear  heterogeneous  "comparison"  material,  with 
effective  properties  that  can  be  characterized  in  terms  of  bounds  and  estimates,  to  obtain 
corresponding  bounds  and  estimates  for  the  effective  properties  of  a  nonlinear  composite.  To  this 
end,  we  introduce  the  quadratic  potential 


6(o.  *)  =  X*w<*>  w«fl  =  7VT  a- +  TTT  a-  ■  (6) 

7*  6m(x)  2xt(x) 

such  that /i(x)  =  ^^(r)(x)/i<r)>0,  and  K,(x)  =  ^£(r)(x)K,('-)>0,  with  the  /2(,)  and  klr) 

r=l  r=l 

constant,  corresponding  to  a  linear  isotropic  composite  with  the  same  phase  distribution  (the  same 
indicator  functions)  as  the  nonlinear  composite. 

Then,  if  the  nonlinearity  in  the  potential  of  the  original  composite  is  stronger  than  quadratic 
as  the  norm  of  the  stress  becomes  large,  it  makes  sense  to  define  the  set  of  functions 

V(r\filr),  klr))  =  max{(7<rS(a)  -  t/(r)(a)},  (7) 

a 

such  that 


V(fi,  k)  =  5>”(x)  kir))  =  max{(7(a,x)  -  C/(a,x)} .  (8) 

r=l  ° 

It  follows  that,  for  all  kir)  >  0  (r  =  1,...,  n)  and  a,  at  each  x  e  Q 

U(a,x)>U(o,x)-V(fi.,k), 

and  hence  that  for  all  /2<r),  klr)  >  0  (r  =  1,...,  n),  and  for  every  a 

U(a)>0(a)-V(fi,k),  (9) 

where 

(7(a)  =  min  (7(a)  (10) 

atS(  5) 

is  the  effective  potential  of  the  linear  composite,  and 

V(fitk)  =  ^clr)V(r\fiir\kir)), 

r»  1 
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expressed  in  terms  of  the  volume  fractions  of  each  phase, 

c{r)  =  \aX(r\x)dV. 

Thus,  if  we  could  compute  0(a)  for  the  linear  composite  in  terms  of  /2(r)  and  k(r),  expression  (9) 

yields  a  family  of  bounds  for  the  effective  potential  of  the  nonlinear  composite,  0(a),  for  every 

choice  of  the  set  of  parameters  j±(r),  k(r)  >  0.  This  family  of  bounds  can  be  optimized  by 
considering 


0_(a)=  max  (0(a)  ~  V(p,k)}. 

Aw.rw>  0 

Then,  evidendy. 


(11) 


0 (a)  >  0_(a) .  (12) 

Usually,  however,  it  is  not  possible  to  find  0(a)  explicidy,  but  instead  bounds  and  estimates  may 
be  available  for  0(a).  If  we  have  a  lower  bound  (such  as  a  Hashin-Shtrikman  lower  bound) 
0_(a)  for  0(a),  such  that 


U(a)>U_(a),  (13) 

then,  replacing  0(a)  by  0_(a)  in  equation  (1 1)  for  0_(a),  still  yields  a  lower  bound  for  0(a); 
alternatively,  an  upper  bound  for  0(a)  is  not  useful  in  terms  of  obtaining  an  upper  bound  for 
0(a).  On  the  other  hand,  if  we  only  have  an  estimate  (such  as  a  self-consistent  estimate)  0t(a) 
for  0(a),  then 


0,(a)=  max  (Ot(a)-V(p,k)}  (14) 

would  provide  only  an  estimate  for  0(a). 

We  note  that  the  prescriptions  (11)  and  (14)  lead  to  convex  expressions  for  the  bounds  and 
estimates  of  the  effective  potential  of  the  nonlinear  composite,  provided  that  the  corresponding 
bounds  and  estimates  for  the  linear  composite  are  convex,  which  is  in  turn  guaranteed  assuming 

that  p  and  k  >  0.  This  is  a  desirable  feature,  because  the  effective  potential  of  the  composite  is 
known  to  be  convex. 

In  reference  [1],  this  general  procedure  was  applied  to  an  isotropic  porous  material  with 
incompressible  behavior  for  the  matrix  with  potential 

U°\a)  =  f(at),  (15) 

where/  is  assumed  to  satisfy  the  condition  that  F(x)  =  f(s)  be  a  convex  function  of  x  -  s2  >  0. 

The  result  of  this  calculation  are  a  Hashin-Shtrikman  (H-S)  lower  bound  for  0 (a),  given  7 

O.(o)  =  c{l)f(s),  (16a) 
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with 


s  =  ^17  V^1  +  %c(2))ct,2  +  %c(2)  a\ , 
and  a  self-consistent  (S-C)  estimate  given  by 

U.(o)  =  cmf(s), 

with 


(16b) 

(17a) 


s  = 


1 

ri-c®/3'| 

c(1) 

[  ‘  4c(1)  "J 

(17b) 


New  Results  for  the  Porous  Rambere-Osgood  Material 

In  this  section,  we  specialize  further  the  results  of  the  previous  section  by  taking  the 
constitutive  behavior  of  the  matrix,  as  given  by  equation  (15),  to  be  described  by  the  relation 

f(s)  =  [f'(s)ds 

where 


f'(s)  =  eA 


/• 

~ 

(  \ 

n 

f  —  y* 

s 

—  + 

s  I 

— 

H(s-oy) 

^oj 

lCTJJ 

\  y  / 

(18) 


Here,  H  is  the  unit  step  function,  n  >  1,  a  is  the  yield  stress,  and  <jo  and  eo  are  used  to 

normalize  the  stresses  and  strains,  respectively,  such  that  the  ratio  oo/eo  corresponds  to  the 
Young’s  modulus  of  the  material.  Then,  the  uniaxial  stress/strain  relation  for  the  matrix  material  is 
linear  up  to  the  yield  level,  and  linear  plus  power- hardening  for  stress  levels  in  excess  of  the  yield 
value.  Note  that  this  function  satisfies  the  convexity  assumption  invoked  above. 

Thus,  the  H-S  bound  and  S-C  estimate  for  the  effective  energy  of  the  porous  material  are 
given  by  expressions  (16)  and  (17),  respectively,  where  /is  specified  by  equation  (18).  From 
these  expressions,  we  can  also  determine,  via  equation  (5),  effective  stress/strain  relations  for  the 
porous  material  of  the  form 


£,  =  |c(1)/'(s)Jk  (19) 

2  do, 

where  s  is  given  by  either  (16b)  for  the  H-S  bound,  or  by  (17b)  for  the  S-C  estimate. 

We  compare  the  new  results  with  the  results  of  reference  [3]  for  the  same  material  in 

Figures  1.  All  the  results  correspond  to  cr^/cro  =  0.2,  n  =  10  and  porosity,  c(2)  =  0.3.  Figure  la 

shows  normalized  plots  of  the  H-S  lower  bounds  and  the  S-C  estimates  for  the  effective  energy  U 

as  functions  of  the  average  effective  stress,  a, ,  for  a  fixed  value  of  the  average  hydrostatic  stress, 

ojo0  =  0.25.  The  dashed  lines  correspond  to  the  new  results  given  by  equations  (16)  and  (17), 
and  the  continuous  lines  correspond  to  the  "best"  results*  of  [3].  The  H-S  lower  bounds  lie  below 


*  The  raw  data  of  [3]  was  not  available,  and  the  plots  of  these  results  in  Figure  1  was  accomplished  by  graphical 
methods.  Thus,  there  may  be  some  small  error  in  the  representation  of  the  results  of  [3]  in  Figures  1. 
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the  S-C  estimates,  as  they  should.  It  is  seen  that,  for  this  case,  corresponding  to  relatively  low 
triaxialities,  the  new  results  are  not  very  different  from  the  results  of  [3],  but  the  new  bound  lies 
slightly  above  the  bound  of  [3],  and  hence  it  is  better.  Figure  lb  shows  similar  normalized  plots 

for  the  case  where  <jJa0  =  l.Q.  In  this  case,  corresponding  to  higher  triaxialities,  the 
improvement  in  the  H-S  lower  bound  is  more  dramatic,  with  the  new  bound  lying  well  above  the 
bound  of  [3].  The  new  S-C  estimate  also  lies  well  above  the  bounds  and  S-C  estimate  of  reference 
[3]  (outside  the  range  depicted  in  the  graph).  This  is  not  surprising,  however,  because  for  such 
large  values  of  n  the  differences  in  results  are  being  exaggerated  by  the  stress  energies  (which  are 

proportional  to  the  stress  to  the  n+1  power!).  Figures  lc  and  Id  show  plots  of  at  versus 
according  to  equation  (19),  corresponding  to  the  bounds  and  estimates  given  in  Figures  la  and  lb, 
respectively.  The  S-C  estimates,  being  less  stiff,  lie  below  the  H-S  bounds.  It  is  interesting  to  note 
that  the  new  results  do  not  show  the  complicated  structure  of  the  results  of  [3]  for  the  larger  value 
of  the  average  hydrostatic  stress.  This  is  because  for  the  new  results,  a  large  enough  average 
hydrostatic  stress  will  saturate  the  linear  response  of  the  material,  in  such  a  way  that  additional 
average  shear  stress  leads  to  a  less  stiff  power-type  response.  By  contrast,  the  stress/strain  relation 
predicted  by  [3]  does  not  seem  to  be  influenced  significantly  in  its  initial  stages  by  the  level  of 
applied  average  hydrostatic  stress,  showing  an  unrealistically  stiff  initial  behavior. 

Figures  2a  and  2b  show  normalized  plots  of  the  H-S  bounds  for  the  effective  energy  U  as 
functions  of  the  average  effective  stress,  at,  in  the  low-triaxiality  range  corresponding  to 


co  = 


«  1,  for  values  of  n  =  3  and  10,  respectively.  The  three  curves  correspond  to  values 


of  c<2)  =  0.25,  0.1  and  0.01,  and  a  value  of  <Jy/<J0  =  0.2.  Higher  porosity  naturally  leads  to  less 
stiff  behavior.  Figures  2c  and  2d  show  normalized  plots  of  the  H-S  bounds  for  the  effective  energy 
U  as  functions  of  the  average  hydrostatic  stress,  on,  in  the  high-triaxiality  range  corresponding  to 
<w»l,  for  values  of  n  =  3  and  10,  respectively.  The  two  curves  correspond  to  values  of 
c(2)  =0.25,  and  0.1. 
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Fig.  2  Normalized  plots  of  (a)  and  (b)  H-S  bounds  for  U  against  G,  in  the  low-triaxiality  range, 

and  of  (c)  &  (d)  H-S  bounds  for  0  against  Gn  in  the  high-triaxiality  range,  for  the  two  values  of 
n  =  3,  10,  respectively.  The  different  curves  correspond  to  three  different  values  of  the  porosity 
c(2)  =  0.25,  0.1  and  0.01,  with  the  higher  porosities  lying  lower  in  the  graph. 
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Abstract 


A  new  variational  method  for  estimating  the  effective  properties  of  nonlinear 
composites  in  terms  of  the  corresponding  properties  of  linear  composites  with 
the  same  microstructural  distributions  of  phases  is  applied  to  an  isotropic, 
incompressible  composite  material  containing  a  brittle  (linear)  and  a  ductile 
(nonlinear)  phase.  More  specifically,  in  this  particular  work  the  prescription  is 
used  to  obtain  bounds  of  the  Hashin-Shtrikman  type  for  the  effective  properties 
of  the  nonlinear  composite  in  terms  of  the  well-known  linear  bounds.  It  can  be 
shown  that  in  some  cases  the  method  leads  to  optimal  bounds. 


Introduction 


PONTE  CASTANEDA  (1990a)  has  proposed  a  new  procedure  for  estimating  the 
effective  properties  of  composite  materials  with  phases  exhibiting  nonlinear 
constitutive  behavior.  The  procedure,  which  is  straightforward  to  implement, 
expresses  the  effective  properties  of  the  nonlinear  composite  in  terms  of  the 
effective  properties  of  a  family  of  linear  composites  with  the  same  distribution 


of  phases  as  the  nonlinear  composite.  Appropriate  references  for  the  linear  theory 
of  composites  are  given  by  the  review  article  of  WILLIS  (1982)  and  by  the 
monograph  of  CHRISTENSEN  (1979).  The  new  procedure  was  applied  in  the 
above  reference  to  materials  containing  a  nonlinear  matrix  either  weakened  by 
voids  or  reinforced  by  rigid  particles.  Estimates  and  rigorous  bounds  were 
obtained  for  the  effective  properties  of  such  materials.  The  Hashin-Shtrikman 
bounds  (obtained  via  the  new  method  from  the  linear  Hashin-Shtrikman  bounds) 
were  found  to  be  an  improvement  over  the  corresponding  bounds  obtained  by 
PONTE  CASTAlsTEDA  and  WILLIS  (1988)  for  the  same  class  of  materials  using  an 
extension  of  the  Hashin-Shtrikman  variational  principle  to  nonlinear  problems 
proposed  by  TALBOT  and  WILLIS  (1985).  Recendy,  WILLIS  (1990)  has  shown 
that  the  Hashin-Shtrikman  bounds  obtained  via  the  new  method  can  also  be 
obtained  by  the  method  of  TALBOT  and  WILLIS  (1985)  with  an  optimal  choice  of 
the  comparison  material.  More  generally,  however,  the  new  procedure  can  make 
use  of  linear  higher-order  bounds  and  estimates  to  yield  corresponding  bounds  and 
estimates  for  nonlinear  materials.  In  fact,  the  new  procedure  can  be  shown  to 
yield  exact  results  for  a  certain  class  of  nonlinear  composites.  This  is  discussed 
in  detail  by  PONTE  CASTANEDA  (1990b). 

In  this  paper  we  apply  the  general  procedure  to  a  composite  containing  a 
brittle  (linear)  and  a  ductile  (nonlinear)  phase.  We  assume  that  the  phases  are 
perfectly  bonded  to  each  other,  incompressible  and  isotropic.  Additionally,  the 
size  of  the  typical  heterogeneity  is  assumed  to  be  small  compared  to  the  size  of 
the  specimen  and  the  scale  of  variation  of  the  applied  loads.  It  is  further  assumed 
that  the  effect  of  the  interfaces  is  negligible,  so  that  the  effective  properties  of 
the  composite  are  essentially  derived  from  the  bulk  behavior  of  the  constituent 
phases.  Both  upper  and  lower  bounds  of  the  Hashin-Shtrikman  type  are  given  for 
the  isotropic  composite  as  functions  of  the  properties  and  volume  fractions  of 
the  phases.  Specific  results  are  given  when  the  behavior  of  the  nonlinear  phase  is 
linear  plus  power-law,  including  the  pure  power-law  case.  Some  of  the  bounds 
are  shown  to  be  optimal  (i.e.,  microstructures  can  be  given  attaining  these 
bounds). 


Effective  Properties 


Consider  a  two-phase  composite  occupying  a  region  of  unit  volume  i2,  with 
each  phase  occupying  a  subregion  £2fr>  (r  =  1,  2),  and  let  the  stress  potential, 

U(a,x),  be  expressed  in  terms  of  the  homogeneous  phase  potentials,  U{r\o), 
via 


f/(o,x)  =  £^',(x){/<°(o), 

r=\ 


(1) 


(2) 


*(,)<x)  = 


1  if  x  e  Cl<r) 

0  otherwise 

is  the  characteristic  function  of  phase  r.  The  phases  are  assumed  to  be 
incompressible  and  isotropic,  so  that  the  potentials  t/(,)(o)  can  be  assumed  to 
depend  on  the  stress  a  only  through  the  effective  stress 


’•-vF*- 


where  S  is  the  deviator  of  o.  Thus,  we  assume  that  there  exist  scalar-valued 
functions  /(,)  such  that 

Uir\a)  =  fir\o.). 

Then,  the  stress  field  <r,  satisfying  the  equilibrium  equations 

aa.j  =  0'  (3) 

is  related  to  the  strain  field  e,  related  to  the  displacement  field  u  via 


1 , 

£V=  2  (“u  + “;.•>* 


(4) 


through  the  constitutive  relation 


V 


du 

do,. 


(5) 


The  commas  in  equations  (3)  and  (4)  denote  differentiation,  and  the  summation 
convention  has  also  been  ut^d  in  equation  (3).  We  assume  that  the  phases  are 
perfectly  bonded,  so  that  the  displacement  is  continuous  across  the  interphase 
boundaries.  However,  the  strains  and,  therefore,  the  stresses  may  be 
discontinuous  across  such  boundaries,  and  hence  equation  (3)  must  be  interpreted 
in  a  weak  sense,  requiring  continuity  of  the  traction  components  of  the  stress 
across  the  interphase  boundaries. 

We  note  that  if  we  let  e  represent  the  rate-of-deformation  tensor  and  u  the 
velocity  field,  the  above  equations  can  be  used  to  model  high-temperature  creep, 
as  well  as  high-rate  viscoplastic  deformations.  Here  we  will  present  our  work  in 
the  context  of  time-independent  plasticity  (deformation  theory),  but  in  view  of 
the  above  comment  the  results  could  be  given  appropriate  interpretations  in 
nonlinear  creep  and  viscoplasticity. 

To  define  the  effective  properties  of  the  heterogeneous  material  we  introduce, 
following  HILL  (1963),  the  uniform  constraint  boundary  condition 


a‘ini  -  a,jnj,  X  e  (9f2,  (6) 

where  dCl  denotes  the  boundary  of  the  composite,  n  is  its  unit  outward  normal, 
and  a  is  a  given  constant  symmetric  tensor.  Then,  the  average  stress  is  precisely 
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ct  =  jo  a(x)dV  ( ') 

and  we  define  the  average  strain  in  a  similar  manner  by 

e  =  [  t(x)dV .  (8) 

JQ 

The  effective  behavior  of  the  composite,  or  the  relation  between  the  average 
stress  and  the  average  strain  then  follows  from  the  principle  of  minimum 
complementary  energy,  which  can  be  stated  in  the  form 

0(a)  =  min  U(a),  (9) 

oeSO) 

where 

U(a)  =  faU(c,x)dV 

is  the  complementary  energy  functional  of  the  problem, 

S(a )  =  {ctI(Tv  >  =  0  in  £2,  and  <Jtjnj  =  on  dO] 
is  the  set  of  statically  admissible  stresses,  and  where  we  have  assumed  convexity 
of  the  nonlinear  potential  £/(a,x).  Thus,  we  have  that 


(10) 


Our  task  will  be  to  determine  bounds  and  estimates  for  U(a),  which,  under  the 
above  assumptions,  is  known  to  be  convex. 


Bounds  and  Estimates 


A  new  variational  principle  for  determining  bounds  and  estimates  for  the 
effective  properties  of  nonlinear  composites  in  terms  of  the  effective  properties 
of  linear  composites  was  proposed  by  PONTE  CASTANEDA  (1990a,b).  In  this 
section,  we  specialize  the  derivation  given  in  PONTE  CASTANEDA  (1990b)  for 
the  case  where  both  phases  are  incompressible,  and  phase  #2  is  linear  so  that 


The  new  variational  principle  is  based  on  a  representation  of  the  potential  of 


the  nonlinear  material  in  terms  of  the  potentials  of  a  family  of  linear 
comparison  materials.  Thus,  for  a  homogeneous  nonlinear  material  with 


"stronger  than  quadratic"  growth  in  its  potential,  t/(o),  we  have  that 
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U(o)>max[U,(o)-V(M)}, 

n>  0 

where 


and 


VOi)  =  ma x{£/.(a)-t/(a)} 


(11) 

(12) 


U.{a)  =  -l-o2.  (13) 

6/X 

is  the  the  potential  of  the  comparison  linear  material. 

To  demonstrate  this  result,  let 

U(a)  =  <t>(s),  (14) 

where  s  =  a ].  Then,  the  Legendre-Fenchel  transform  of  the  scalar- valued 
function  0  is  given  by 

<p‘(a)  =  max{ar  -  0(s)},  (15) 

j>0 

where  a  is  assumed  to  be  positive.  A  well-known  result  in  convex  analysis 
(VAN  TIEL  1984,  §6.3)  is  that 

0(s)>max{.ra-0*(a)},  (16) 

a>0 

with  equality  if  0  is  a  convex  function  of  its  argument.  With  the  identifications 
s  =  and  a  =  (6/t)~\  we  can  see  that  (11)  and  (12)  are  but  simple  re¬ 
statements  of  (16)  and  (15),  respectively.  In  particular, 


(17) 


To  derive  the  new  variational  principle,  we  apply  (1 1)  to  the  nonlinear  phase 
#1,  and  make  use  of  the  result  in  the  complementary  energy  principle  (9).  Thus, 
after  some  manipulations,  we  Find  that 


t/(5) 


>  max  \u 


(a)-  \Vwm"')dV 


}• 


(18) 


where 


t/„(a)=  min  t/0(a), 

0€$(5) 


(19) 


U0(a,x)  =  ^x(,)^Wlr\o), 

r-l 


^,)(o)=7J<o^- 


6n 


and 
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Note  that  the  comparison  linear  material  agrees  with  the  actual  material  in  phase 
#  2  (which  is  linear).  In  the  above  derivation,  we  note  that  the  comparison 

moduli  /z(1)  are  functions  of  x,  since  the  stress  field  a  will  also  in  general  be  a 

function  of  x  within  phase  #1.  If  we  assume  that  Uw(a)  is  "strongly  convex” 
(i.e.  if  <(>  is  convex),  then  we  have  equality  in  (11),  and  hence,  usually,  equality 
in  (18).  However,  if  the  conditions  for  equality  are  not  met,  relation  (18)  still 

provides  a  useful  lower  bound  for  0(a).  An  detailed  derivation  of  this  result, 
discussing  the  precise  conditions  for  equality,  is  given  in  PONTE  CASTANEDA 
(1990b). 

The  variational  principle  described  by  (18)  roughly  corresponds  to  solving  a 
completely  linear  problem  for  a  heterogeneous  material  with  arbitrary  moduli 
variation  within  the  nonlinear  phase,  and  then  optimizing  with  respect  to  the 
variations  in  moduli  within  the  nonlinear  phase.  Thus,  one  can  think  of  the 
nonlinear  material  as  a  "linear"  material  with  variable  moduli  that  are  determined 
by  prescription  (18)  in  such  a  way  that  its  properties  agree  with  those  of  the 
nonlinear  material. 

This  suggests  that  if  the  fields  happen  to  be  constant  over  the  nonlinear 
phase,  then  the  variable  moduli  M°’(x)  can  be  replaced  by  constant  moduli  nw. 
More  generally,  however,  we  have  the  following  lower  bound  for  U(a) 


0  (a)  =  maxft/^o)  -  c(,V<,>(/i0>)) .  (20) 

M<1’>  0 

where  cw  is  the  volume  fraction  of  phase  #1.  The  result  in  this  form  is  a  special 
case  of  a  more  general  result  first  derived  by  PONTE  CASTANEDA  (1990a),  when 
only  one  of  the  phases  is  nonlinear,  and  the  other  one  is  linear. 

We  note  that  the  prescriptions  (18)  and  (20)  lead  to  convex  expressions  for 
the  bounds  and  estimates  of  the  effective  potential  of  the  nonlinear  composite, 
provided  that  the  corresponding  bounds  and  estimates  for  the  linear  composite  are 
convex.  This  is  a  desirable  feature,  because  the  effective  potential  of  the 
composite  is  known  to  be  convex. 


Application  to  Hashin-Shtrikman  Bounds 


HASHIN  and  SHTRIKMAN  (1962)  prescribed  bounds  for  the  effective  moduli  of 
linear-elastic,  isotropic  composites,  depending  only  on  the  volume  fractions  of 
the  phases.  When  there  are  only  two  phases,  these  bounds  have  been  shown  to 
be  optimal  (i.e.,  microstructures  can  be  given  that  simultaneously  attain  the 
bounds  for  the  shear  and  bulk  modulus)  by  FRANCFORT  and  MURAT  (1987). 
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8  «  e 

Figure  1.  Rank-2  laminate. 


Their  construction  made  use  of  iterated  laminates  for  which  the  effective 
properties  can  be  computed  exacdy.  Such  materials  are  obtained  by  layering  the 
two  constituent  phases  to  obtain  a  rank-1  laminate;  the  resulting  material  is  once 
again  layered  (in  an  arbitrary  direction)  with  one  of  the  original  phases  in  a 
smaller  lengthscale.  This  procedure  can  obviously  be  iterated  n  times  to  obtain  a 
rank-n  laminate.  In  general  such  materials  will  be  anisotropic,  but  by  choosing 
appropriately  the  layer  orientations  at  the  different  layering  operations,  it  is 
possible  to  obtain  an  isotropic  composite,  and  its  properties  coincide  with  one  of 
the  Hashin-Shtrikman  (H-S)  bounds  depending  on  which  constituent  phase  is 
selected  to  play  the  role  of  the  matrix  material.  Figure  1  depicts  a  rank-2 
laminate  (not  to  scale)  with  phase  #2  as  the  matrix  phase. 

For  the  special  case  of  incompressible  materials,  when  there  is  only  one 
modulus  for  the  composite,  the  H-S  upper  bound  for  the  effective  shear  modulus 
can  be  expressed  in  the  form 


where 


/r ' 

nm 


if  nw>nm 
if  <  fim 


a(M°VJ))  = 


2c(V2,+(3  +  2  c<JV!> 
(2  +  3c(2><2,  +  3c<V> 


(21) 


(22) 
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and 


P(n"\nm)= 


2cmnll)+{3+2  cm)nw 
(2  +  3c<!>cl)  +  3c(V2>  ' 


(23) 


The  corresponding  H-S  lower  bound  is  obtained  by  interchanging  the  expressions 
in  (21)  for  the  upper  bound  (and  keeping  the  conditions  on  the  shear  moduli  fim 
and  nw  fixed). 

The  above  H-S  upper  bound  for  the  effective  shear  modulus  yields  a  lower 
bound  for  the  potential  of  the  linear  material  U„ .  This  information  can  be  used  in 
combination  with  prescription  (20)  to  yield  a  H-S  lower  bound  for  the  potential 


of  the  nonlinear  material  U.  On  the  other  hand,  upper  bounds  for  U0  do  not 


necessarily  generate  upper  bounds  for  0 . 

The  result  for  the  lower  bound  on  0  depends  on  which  of  the  two  branches 
of  (21)  is  used  in  conjunction  with  (20).  If  nm  >  /t<2>.  then  the  average  effective 
stress  a,  must  be  such  that  the  condition 


3 mmf\o,)<G„  (24) 

is  satisfied  (usually  when  the  average  shear  stress  is  small  enough).  Here,  for 
simplicity,  we  have  made  the  identification  f0)  =  f .  The  corresponding  form  of 
the  bound  is  then 

0.(3)  =  /,(?.),  (25) 

where 


/,(<t.)  =  c(1,/(5)+m!J> 

and  s  solves  the  equation 


j>3c-)  fa, 


(ns)f 


(26) 


cH)  +/r<2,(2  +  3c<2>)^-^  =  — 


(2  +  3c(2>)(s/g.)  -2 
3cl2> 


(27) 


On  the  other  hand,  if  /x(I>  <  nw,  then  the  average  effective  stress  a,  must  be 
such  that  the  condition 

3/r<2>/'(g.)  >  a„  (28) 

is  satisfied  (i.e.,  when  the  average  shear  stress  is  large  enough),  and 

0(o)  =  f2(a.), 


(29) 
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where 


(3  +  2ca))^  +  cm(2  +  3c{1Y  -  \Ocwsa. 

A(*.)  =  cwf(s)  +  ~ - 1  1 - 

and  s  solves  the  equation 


9cm(J.mf'{s)  =  5 a,  -  (2  +  3ca))s . 
The  corresponding  stress/strain  relations  have  the  form 

E=|m)s, 

where 


(30) 

(31) 

(32) 


2  (3  +  2c(1))g,-5c(1>s 
Mt>  3  6  cmnm 


(33) 


but  filer.)  does  not  have  a  simple  expression. 

In  general,  we  do  not  expect  the  above  lower  bounds  for  0  to  be  optimal.  In 
fact,  expression  (25)  does  not  yield  an  optimal  bound  if  condition  (24)  is 
satisfied.  However,  it  is  shown  in  PONTE  CASTANEDA  (1990b)  that  if  condition 
(28)  is  satisfied,  then  the  bound  (29)  is  optimal.  This  is  because  the  same 
microstructure  attaining  the  linear  bounds  can  be  also  shown  to  attain  the 
nonlinear  bound;  the  reason  being  that  the  fields  are  constant  in  the  (nonlinear) 
inclusion  phase,  and  hence  expressions  (20)  and  (18)  are  identical.  Similar 
observations  have  been  made  by  KOHN  (1990)  in  a  similar  context  (starting  from 
the  Talbot-Willis  nonlinear  variational  principle)  and,  independently,  by  PONTE 
CASTANEDA  ( 1990c)  i  the  context  of  conductivity. 

Conversely,  in  general,  we  do  not  expect  that  interchanging  conditions  (24) 
and  (28)  would  turn  expression  (25)  and  (29)  into  upper  bounds  for  the  nonlinear 
potential  0.  This  is  contrary  to  the  corresponding  operation  for  the  linear 
composite.  All  that  can  be  said,  however,  is  that  expression  (29)  is  an  estimate 
for  the  upper  bound  for  0  if  condition  (24)  is  satisfied  and  that  expression  (25) 
is  an  estimate  for  the  upper  bound  for  0  if  condition  (28)  is  satisfied.  Both  of 
these  estimates  are  expected  to  get  progressively  better  with  weaker 
nonlinearities. 


Application  to  Power-Law  Behavior 


In  this  section,  we  specialize  further  the  calculations  of  the  previous  section  by 
taking  the  constitutive  behavior  of  the  nonlinear  phase  to  be  governed  by  a  linear 
plus  power  relation 
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(34) 


Note  that  the  case  n  ->  «  corresponds  to  pure  power-law  behavior,  and  the 
limits  n  — » 1  (in  addition  to  /t  — ►  °°)  or  rj  — » <»  correspond  to  linear  behavior. 

The  conditions  (24)  and  (28)  determining  the  appropriate  branch  of  the 
bound  specialize  to 


(35) 


and  the  opposite  inequality,  respectively.  The  first  condition  guaranteeing  that 
(25)  is  a  lower  bound  (and  correspondingly  that  (29)  is  an  estimate  for  the  upper 
bound)  corresponds  to  small  enough  average  stress  on  the  composite. 
Alternatively,  the  second  condition  (with  >  instead  of  <)  corresponds  to 
sufficiently  large  average  stress.  Note  that,  if  ni2)/n  >  1 ,  condition  (35)  can 
never  be  satisfied  and,  conversely,  the  alternative  condition  is  always  satisfied. 
This  condition  ensures  that  the  difference  between  the  potential  of  phase  #1  and 
that  of  phase  #2  is  convex.  Here,  we  will  consider  two  cases:  one  case,  meeting 
this  condition,  with  /fi  =  2,  and  the  other  with  IfJ-  =  0,  corresponding  to 
the  pure  power-law  case. 

The  results  for  the  bounds  (25)  and  (29)  specialized  to  the  case  when  (34) 
holds  can  be  expressed  in  the  form: 


0(o) 

Um(a) 


(36) 


where  the  precise  form  of  F  depends  on  whether  (25)  or  (29)  applies,  and 


plays  the  role  of  the  independent  variable,  with  cm  and  n, 

serving  as  parameters. 

Results  for  the  upper  and  lower  bounds  for  0  are  given  in  Figures  2  and  3 
for  the  case  where  =  0,  and  in  Figure  4  for  the  case  where  fim/n  =  2.  In 
the  first  case,  condition  (35)  determining  whether  (29)  is  an  estimate  for  the 
upper  bound  or  an  optimal  lower  bound,  and  whether  (25)  is  an  estimate  for  the 
upper  bound,  or  a  non-optimal  lower  bound,  simply  reduces  to  the  condition  of 


whether  the  independent  variable  '  is  less  or  greater  than  unity.  For 

that  reason,  we  give  results  emphasizing  the  small  stress  and  large  stress 
domains,  separately,  in  Figures  2  and  3,  respectively. 

In  each  plot  we  have  three  sets  of  curves  corresponding  to  three  values  of 
c<2>  (0.1,  0.5  and  0.9).  Additionally,  we  show  the  limiting  cases  corresponding 
to  cm=  0  and  cm=  1.  These  limiting  curves  appear  as  straight  lines,  one  with 
with  variable  slope  depending  on  the  value  of  n  and  nm/n,  and  the  other  with 


U/U<2)  1 


0.8 

0.6 

0.4 

0.2 
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(n<2>/n(1))jen'1 

Figure  2(a).  Plots  of  the  bounds  for  the  effective  energy  of  the 
composite  as  functions  of  the  average  stress  (appropriately  normalized) 

for  =  0  and  n  =  3  (small  stress). 


zero  slope  (value  equal  to  unity),  respectively.  The  intermediate  sets  of  curves 
correspond  to  the  upper  and  lower  bounds. 

In  Figure  2,  depicting  results  for  two  values  of  the  nonlinearity  parameter 
(n  =  3  and  10),  the  continuous  line  corresponds  to  the  estimate  for  the  upper 

bound  (for  0),  and  the  dashed  line  corresponds  to  the  rigorous  lower  bound.  In 
Figure  3,  showing  also  results  for  the  same  two  values  of  the  nonlinearity 
parameter,  the  continuous  line  corresponds  to  the  optimal  lower  bound,  and  the 
dashed  line  is  an  estimate  for  the  upper  bound.  For  this  value  of  M<2)/m.  the 
upper  and  lower  bound  coalesce  when  the  value  of  the  independent  variable 
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(ii(2)/n(1)>o(.n'1 

Figure  2(b).  Plots  of  the  bounds  for  the  effective  energy  of  the 
composite  as  functions  of  the  average  stress  (appropriately  normalized) 

for  n"’ /n  =  0  and  n  =  10  (small  stress). 

[fim approaches  unity.  In  the  linear  case,  this  behavior  corresponds  to 
the  limit  of rr.rirM  cf  the  pha'es  appro-,*. hint*  f>ach  other.  More  generally, 
assuming  that  /xa>//*  is  less  than  unity,  there  is  a  value  of  the  independent 
variable  (i.e.,  an  average  stress  level)  at  which  the  bounds  are  equal,  and  hence 
the  effective  energy  of  the  composite  is  known  exactly.  This  phenomenon  is 
related  to  the  lack  of  convexity  of  the  difference  between  the  potentials  of  the 
nonlinear  and  linear  phases. 

In  Figure  4,  depicting  results  for  the  same  two  values  of  the  nonlinearity 
parameter,  the  continuous  line  corresponds  to  the  optimal  lower  bound  (for  0), 
and  the  dashed  line  corresponds  to  the  estimate  for  the  upper  bound.  In  this  case. 
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Figure  3(a).  Plots  of  the  bounds  for  the  effective  energy  of  the 
composite  as  functions  of  the  average  stress  (appropriately 

normalized)  for  n'"/n  =  0  and  n  =  3  (large  stress) 


with  a  convex  difference  between  the  nonlinear  and  linear  potentials,  there  is  no 
value  of  the  independent  variable  for  which  the  upper  and  lower  bound  are  equal. 

Both  in  Figures  3  and  4,  we  observe  that  the  lower  bound  approaches  a 
straight  line  with  zero  slope  and  the  upper  bound  approaches  a  straight  line  with 
slope  depending  on  the  value  of  n  (smaller  for  larger  n).  This  is  consistent  with 
the  following  asymptotic  behaviors  for  the  lower  and  upper  bounds 

jm.Oixp,  (37) 

Um(a)  cm 


x 
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Figure  3(b).  Plots  of  the  bounds  for  the  effective  energy  of  the 
composite  as  functions  of  the  average  stress  (appropriately 

normalized)  for  =  0  and  n  =  10  (large  stress) 
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Figure  4(a).  Plots  of  the  bounds  for  the  effective  energy  of  the 
composite  as  functions  of  the  average  stress  (appropriately 

normalized)  for  11''  j n  =  2  and  n  =  3. 

respectively.  These  two  behaviors  correspond  physically  to  the  cases  of  a  linear 
matrix  with  voids  and  a  power-law  matrix  with  rigid  inclusions  (studied  by 
PONTE  CASTANEDA,  1990a),  respectively.  The  reason  for  these  behaviors  is  that 
the  lower  bound  (for  0 )  corresponds  to  putting  the  stiffer  material  in  the  matrix 
phase  and  the  less  stiff  material  in  the  inclusion  phase  (and  viceversa  for  the 
upper  bound).  Clearly,  for  large  enough  stresses,  the  linear  phase  is  stiffer  than 
the  nonlinear  phase. 

We  note  that  accurate  numerical  calculations  of  the  potential  of  a  power-law 
matrix  with  spherical  rigid  inclusion  have  yielded  results  of  the  form  (38)  with 

(0  =  1  +  jc<2)  asc(2)-»0. 
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Figure  4(b).  Plots  of  the  bounds  for  the  effective  energy  of  the 
composite  as  functions  of  the  average  stress  (appropriately 

normalized)  for  /t“’  j n  =  2  and  n  =  10. 


where  g(n)  is  such  that  g(l)  =  5/2,  g(3)  »  3.21  and  g(10)  =  6.09  (LEE  and  MEAR, 
1990),  and  g(n)  ->  0.38 n  asn  (HUTCHINSON,  1990).  These  results  do  not 
compare  very  favorably  with  the  corresponding  results  from  (39):  5/2, 4.00,9.25 
and  0.75 n,  but  it  should  be  recalled  that  these  results  correspond  to  the  case  for 
which  we  do  not  have  a  rigorous  bound  (it  is  simply  an  estimate  of  the  bound). 
None  the  less,  the  results  of  (38)  with  (39)  may  provide  reasonable  estimates  for 
larger  values  of  the  volume  fraction  of  the  linear  phase. 
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